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PREFACE 

This monograph is intended to provide a critical review and 
analysis of recent studies of the electromagnetic structure of 
nucleons. The first chapter is devoted to a precise definition of 
what is meant by the electromagnetic form factors. The 
different types of experiments which have given information 
on these form factors are discussed in Chapter II with 
emphasis on the theoretical assumptions and approximations 
made in analysing the experimental numbers in terms of 
form factors. Chapter III is devoted to the recent dispersion 
theoretic analyses of the form factors. The theoretical 
development here has the two goals of serving as a first 
introduction to the dispersion methods for physicists at the 
graduate level and of indicating clearly the various approxima 
tions made at present in implementing the dispersion theory 
approach. In Chapter IV the validity of quantum electro 
dynamics is discussed and it is shown how possible deviations 
in this theory alter the nucleon form factor analyses. Two 
appendices are included with formal developments. 

We wish to thank the Air Force Office of Scientific Research, 
of the United States Air Force, for partial support in prepara 
tion of the manuscript. 

S. D. D. 
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CHAPTER I 

DEFINITION OF ELECTROMAGNETIC 
STRUCTURE OF NUCLEONS 

THE first step in any discussion of tlie electromagnetic 
structure of the nucleon is of necessity a precise definition of 
what is meant by such structure. For purposes of introducing 
and illustrating this definition, it will be convenient to fix 
attention, for the moment, on the scattering of electrons from 
protons; it is this experiment, after all, which has produced 
the present information on the proton's electromagnetic 
structure. 

To lowest order in e, the electric charge, a single virtual 




FIG. I.I. Feynman diagram for electron-proton scattering to order e 2 . 
The proton is treated as a point Dirac particle. The line represents 
a proton, the heavy line, an electron, and the wavy line, a photon. 

photon is exchanged between the electron and the proton. The 
Feynman diagram representing this process is illustrated in 
Fig. I.I. 
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The transition amplitude corresponding to this diagram may 
be written [1] 



(1.1) 

where p fi p i and p ( /\ pW are the final and initial 4-momenta 
of the proton and electron, respectively, q^ = (PfPi)^ = 
(Pi^-Pf^)f, is "kh- 6 4-momentum transferred from electron to 
proton, and jjf\x), ffi(%) are the electromagnetic current 
operators of the proton and electron. The factor q~ 2 is the 
propagator for the virtual photon which is exchanged [2]. 
According to quantum electrodynamics, the current operator 
for the electron is given by, 

j ( $(x) = e^ e (x)y^ e (x) 

where y) e (x) is the electron field, and e is the unrenormalized 
charge. Thus, to lowest order, 



= V( ( 

Similarly, according to pure quantum electrodynamics, the 
proton current is 

j<f\x) = e V(*)y^(*) (1.3) 

so that, again to lowest order in e , 



The lowest order Feynman amplitude is then obtained by 
combining equations (1.1), (1.2), and (1.4). Now the basic 
assumption involved in deducing equation (1.4) is that the 
coupling of the proton to the virtual photon is described 
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simply by quantum electrodynamics, in which the proton is a 
point particle of unit charge and the usual Dirac magnetic 
moment. It is clear, however, that this assumption is not at 
all valid in particular the proton couples strongly to mesons 
and other particles, so that equation (1.3) for the electro- 
dynamic current of the proton is certainly not correct. For 
example, since there is presumably a cloud of virtual mesons 
around the proton, one should add the meson current to j^, 
i.e. one should write [3] 



where y(x) is the meson field. Similarly more terms should 
be added corresponding to other particles which can appear 
virtually. All these virtual particles, then, produce effects 
which change the structure of the vertex at which the photon 
is absorbed by the proton; that is they change the form of 
equation (1.4). That such effects are not negligible is illus 
trated, for example, by observing that the proton's magnetic 
moment, 2.79e/2JSf, differs greatly from the Dirac moment, 
e/2Jf. 

One must therefore expect deviations from the result of 
equation (1.4). These deviations will, in general, be of two 
kinds. Erst, the fact that the proton's moment is not equal 
to the Dirac moment means that there will be a Pauli term of 
the form a^F ( $ in addition to the simple >yl e) term already 
included [4]. Second, because of the cloud of virtual mesons 
about the proton, the coupling of the photon to the proton may 
take place over an extended region in space thus one should 
replace 



in which the photon is absorbed on a point proton, by 
something like 



d V 
This may be conveniently rewritten as 
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where 



. (1.5) 

Physically, one might interpret this by saying that if the spatial 
extension of the proton is, for example, due to a virtual meson 
cloud, then the photon could be absorbed at x l by a virtual 
meson emitted from the proton at x. F(x x') might therefore 
be expected to have the 'size' of the meson cloud around the 
proton. If these alterations are incorporated into equation 
(1.4), it becomes 




(1.6) 

Here a different spatial extent is associated with the original 
convection current coupling y^J?*, and with the subsequently 
introduced Pauli term a^F^. Note that the charge e Q has 
been included in F^. 

This form has been obtained on the basis of physical ideas 
about what effects could be produced by the presence of a 
meson cloud around the proton. It is in addition, however, the 
most general form which is allowed for the coupling of a 
photon to a physical proton by the requirements of gauge 
invariance and Lorentz invariance. This can be seen as 
follows. Consider the matrix element 



The translational invariance of the theorists' world implies 



where P is the total energy momentum operator; thus 
equation (1.7) becomes 
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Now the matrix element (^^(O)!^) must be of the form 

where 0^ is a Lorentz 4-vector, and is also a 4 x 4 matrix in the 
spinor space of the proton. The factor ^/(^EfE^ appears here, 
incidentally, only for normalization purposes. 

The matrix element is a function only of the 4-vectors p f 
and p i9 which are subject to the restrictions pj = pf = If 2 ; 
from these one can construct only one scalar variable, which 
is chosen to be (p f p t )* = # 2 . The operator 0^ then, must be 
constructed from the scalar g 2 , the vectors p fftt and p f 9 or 
equivalent^ the vectors q^ and P^ = (p f +pi\, and the Dirac 
y matrices. Since 0^ is a 4 x 4 matrix, it can always be written 
as a linear contribution of the 16 matrices 



since, in addition, it is not a pseudovector, the two y 5 terms 
cannot appear. 

The coefficients of the remaining terms are functions of 
whatever scalars can be formed; that is, of q 2 , jf f and 3^. By 
use of the commutation rules of the Dirac matrices, $ f and 3^- 
can be moved to the left and right of the expression, respec 
tively. But & 9 fl f = u Pf M 9 and jf^ = Mu 9 ; 9 thus the tf f 
and tft scalars can be eliminated and the most general form for 

Now the following identities are easily verified: 

(u n \ia rin P^u m .) = f!L I a..\u~.)* 



for free physical proton states which satisfy the Dirac 
equation ($M)u p == 0. Hence 0^ can be reduced to the form 



0, = a'( 



where a', c 7 , and rf' are certain linear combinations of a, 6, c, cZ, 
and e. 
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Up to this point only Lorentz invariance and the fact that 
the initial and final protons satisfy the free particle Dirac 
equation have been used. It is still necessary to impose the 
requirement of gauge invariance. This implies the continuity 
equation Vj^ = 0. Translating this into a condition on 0^ 



therefore 
= 




The last term here is directly zero; the second vanishes 
because 



Thus gauge invariance requires 

gV(# 2 ) = 0. 

Since in general, for virtual photons, g 2 ^ 0, it follows that 
a'(g 2 ) = 0. It is possible to show that a'(0) = as well, by 
observing that the matrix element must be invariant under the 
transformation of interchanging p f and p i9 and reversing the 
order of all y matrices. This is most easily seen by an inspection 
of the relevant Feynman diagrams. The terms y^ and a^ v q v 
are invariant to this, but the term q^ changes sign. 

The conclusion is that the most general form consistent with 
Lorentz invariance and gauge invariance, and with the fact 
that the two protons p f and p i satisfy the free Dirac equation, 
is 



with the scalars c' and d f denoted here as F^ and iF 2 . Purely 
from these invariance arguments, then, it is clear that the 
form given in equation (1.6), on the basis of physical arguments 
is, in fact, the most general one possible. It is easy to see, in 
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addition, from the fact that j^ is a Hermitiaii operator, that 
the J"s are real: explicitly, 



where g' = #< p, = g, so that g' 2 = g 2 . 




FIG. 1.2. Diagram for electron-proton scattering to order e*. The 
blob includes effects of all strongly interacting virtual particles. 

The electron-proton scattering amplitude, to lowest order in 
the electric charge, is then given by 



x 

This includes all possible effects due to clouds of virtual 
particles around the proton, and in fact, due to any contribu 
tion to the proton vertex which does not violate Lorentz 
invariance and gauge invariance. 

The amplitude (1.9) may be interpreted as corresponding to 
a diagram like Fig. 1.2 instead of Fig. 1. 1 . The blob around the 
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proton-photon vertex is supposed to indicate that no attempt 
is made to specify in detail how the interaction occurs that is, 
the blob is meant to summarize all effects of any number of 
virtual particles as indicated in Pig. 1.3. 




-t- 



ywv, 




FIG. 1.3. Examples of diagrams which contribute to the blob of 
Fig. 1.2. The dashed lines represent pions. 

The cross-section resulting from the amplitude of equation 
(1.9) is the Rosenbluth cross-section [5], 



da 



cos 2 (|(9) 



sn 



where 



(2E sin 



0177 



(1.10) 



and U denotes the incident energy and the scattering angle 
of the electron in the laboratory system. 
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If this formula is used to interpret the experimental results 
for electron-proton scattering, then F^ and F 2 may be 
measured as functions of the invariant momentum transfer of 
the scattering. By definition, F^O) = e, the renormalized 
physically observed charge, and F 2 (0) = ^ = g(e/2M) = 
1-79/jjg, the anomalous magnetic moment [6], 

A point particle of charge e and total magnetic moment 
(ff + l)/*B is a particle for which F^q*) = e and ^ 2 (# 2 ) = ^^ 
for all values of g 2 ; this is easily seen from equation (1.5). Thus 
by definition, a particle has electromagnetic structure i.e. is 
not a point particle if and only if the functions F^(q 2 ) and/or 
F 2 (q 2 ) are not constant. This is the precise definition of 
electromagnetic structure. The functions F^ and F 2 are called 
respectively, the charge and moment form factors of the 
proton. 

It should be noted that the definition of F t and F 2 in 
equation (1.8), is quite independent of the particular process 
of electron-proton scattering, but involves only the interaction 
of a real proton with a virtual photon. The same vertex 
appears in many other processes as well; it was only to make 
the physical meaning of the form factors a little clearer that 
their definition was introduced through the discussion of 
electron-proton scattering. 

The functions F^q 2 ) and F 2 (q 2 ) as defined in equation (1.8) 
from consideration of electron proton scattering, are defined 
only for negative i.e. spacelike values of g 2 , because q* = 
(Pf-Pi) 2 = SJfa-^p, .^<0. The definitions of F^ and F 2 
may be extended to include positive i.e. timelike values of 
q 2 by writing 



where <pp /( - } |j M (0)|0> is the amplitude to create a proton- 
antiproton pair from a virtual photon. Here $>' is the proton 
4-momentum, p is the antiproton 4-mornentum, q = p' +p, 
and (pp f( ~ } \ is a proton-antiproton state with ingoing boundary 
conditions [7]. This definition gives F^ and F 2 for 
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It is possible to show that this definition of F- L and F 2 in the 
region q z ^4=M 2 coincides with the analytic continuation of the 
form factors defined in equation (1.8), treated as functions of 
the complex variable g 2 . It will turn out that it is for many 
purposes more convenient to think of J\ and F 2 in terms of 
this pair production process. 

For spacelike q, where F^ and F % may be thought of as 
related to the absorption of a virtual photon on a proton, it is 
possible to choose a co-ordinate system in which q has no time 
component. This is the barycentric system in which the 
incident (or outgoing) electron and proton have equal and 
opposite momenta. In this system, it is conventional to define 
spatial charge and moment densities by 



i(-tf)e-*. (1.11) 

It should be emphasized that this definition depends on a 
particular choice of Lorentz frame, in which the proton is not 
stationary, and therefore the relation of these densities to any 
real physical extension of the proton is quite unclear. 
Note that 



The mean square radii of the spatial distributions are defined 

by 

<r> lfi = JV6 lt ,(r) 



Therefore, for small q z , one gets, in the special co-ordinate 
system: 



This may clearly also be written covariantly in any co 
ordinate system: thus 
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The mean square radius, then, just measures the slope of the 
form factor as q 2 -> 0; if the form factor is a very smooth 
function of q 2 , this slope may give a way of simply summarizing 
its properties. If, however, the form factor does not have a 
fairly constant slope, the value of the rms radius is quite 
independent of its overall behaviour and therefore may in 
general not be a very useful concept. 

Up to this point the discussion has been limited to the form 
factors and electromagnetic structure of the proton. It must 
now be expanded to recognize the existence of neutrons as well. 
The first step in accomplishing this is to affix an index P, for 
proton, onto all the form factors, densities, and radii so far 
defined. 

The same considerations given above for the proton may 
also be made for the neutron; explicitly, define 



where q = n f n iy and n f , n i are the 4-momenta of the neutron 
after and before absorbing the virtual photon. Equivalently, 
for timelike g a >4Jf 2 , the pair production amplitude is 



where q = n' +n, and n' 9 n are the 4-momenta of the neutron 
and antineutron produced by the photon. Then the static 
limits for the neutron fix [8] 

jFf (0) = 



These equations, together with the analogous ones for 
protons, may be combined into one as follows: write 



)cr /lA K) (1.11) 

where p t and p t now refer to the 4-momenta of the final and 
initial nucleons, r s is the 2-component of the isotopic spin 
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operator for the nucleon, and has the value + 1 if p f , p t refer to 
a proton, and 1 if p f , p i represent a neutron. F[ and F% are 
called the isotopic scalar charge and moment form factors, 
Fl and F\ are called the isotopic vector charge and moment 
form factors. 

Evidently the following correspondence obtains: 



Therefore 

Fl(Q) = F*(Q) = e/2 



= Fi>+Fl 

= FI-FI 

= jP| J|. (1.12) 



This decomposition into isotopic scalar and vector form 
factors is useful in the theoretical discussions which are based 
on the assumption of charge independence in the strong 
couplings between mesons and nucleons [7]. In all theoretical 
analyses presented in the succeeding chapters, charge inde 
pendence, i.e. conservation of isotopic spin, will be assumed. 



CHAPTER II 

THEOKETICAL APPROXIMATIONS AND 

ASSUMPTIONS IN OBTAINING NUCLEON 

FORM FACTORS FROM EXPERIMENT 

It was remarked in Chapter I that the definition of the form 
factors FI and F 2 involves only the interaction of a real 
nucleon with a virtual photon. Thus any experimental 
arrangement which replaces the 'question box' of Fig. II. 1 
by a known interaction probes electromagnetic structure of 
the nucleon. Five types of experiments which have been 
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carried out to determine the nucleon form factors are discussed 
in this chapter. These are: 

1. Elastic electron-proton scattering for information on 
proton structure. 

e+P -+e + P 

2. Neutron-electron interaction for information on neutron 
structure. 

N+e -> N+e 




FIG, II. 1. General process for measuring nucleon form factors. 

3. Inelastic electron deuteron scattering for information on 
neutron structure. 

e+D -+e'+N+P 

4. Elastic electron-deuteron scattering for information on 
neutron structure. 

e+D ->e+D 

5. Electron production of pions from protons primarily for 
information on neutron structure. 



e+P 



Theoretical assumptions and approximations which are re 
quired in order to extract values for F^q*) and jF 2 (2 2 ) from 
these experiments are analysed in what follows. 
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1 [9]. For elastic electron-proton scattering the 'question 
box' of Fig. II. 1 is replaced by an electron vertex, Fig. 1.2. 
Hence this process gives the form factors of the proton, jPf 
and F%, directly so long as it is adequate to treat the electro 
magnetic interaction between electron and proton in lowest- 
order perturbation theory, i.e. with one-quantum exchange. 
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FIG. II. 2. Curves for the independent determination of the charge 
and moment form factors of the proton. (Taken from B. Hofstadter; 

see [9].) 

If the validity both of the Born approximation treatment of 
the electromagnetic interaction and of quantum electro 
dynamics is assumed, then the experiments on elastic 
scattering of electrons from protons must give a cross-section 
of the form of the Eosenbluth formula. Comparison of this 
formula with the experiments at various energies and angles 
then yields F^(q 2 ) and F^(q 2 } directly for space-like momentum 
transfers. Evidently two measurements at different energies 
and angles, but fixed q 2 , provide separate information on F l 
and F 2 . 

Charge scattering, or F 17 is of greater importance at small 
angles and moment scattering at large angles. Fig. II. 2 
shows the sensitivity with which two measurements at 
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different energies E and angles 6, corresponding to the same q 2 , 
determine F l and F 2 independently. Drawn in Fig. II. 2 are 
portions of ellipses in the F-F% plane which correspond to a 
given cross-section value for fixed E and 6. A single cross- 
section measurement gives possible values of ( F v J? 2 ) along 
the arc of the ellipse corresponding to a given value of E and 6. 
The intersection with a second ellipse [10] for different values 
of E and 6 but for the same # 2 then uniquely specifies ( F I} F%). 
The figure is drawn for a particular pair and the family of 
three ellipses shows the latitude resulting from an experi 
mental uncertainty of 10 per cent. If one continues this 
procedure and performs a third (or more) measurement with 
yet another pair of energy and angle values with the same <? 2 
there results a consistency check on the analysis. 

Since the form factor analysis of the electron-proton 
scattering experiments are based on the Bosenbluth formula, 
it is appropriate to analyse next its quantitative validity, 
within the framework of quantum electrodynamics. The only 
possible corrections are then due to higher-order electro 
magnetic effects. 

First of all there is the radiative correction to the electron 
current as calculated by Schwinger [11]. The Feynman 
diagrams corresponding to this correction appear in Fig. II. 3. 
The emission of soft photons in the last two diagrams is allowed 
as consistent with the experimental energy resolution. For 
the high resolution used in the work of Hofstadter et al. [9], 
this correction may be as large as 15%~20% and it is taken 
into account in reducing the experimental data for comparison 
with the Eosenbluth formula. The analogous radiative 
correction to the proton current is negligible in present 
experiments due to the large proton mass. 

A second class of possible corrections to the form factor 
interpretation corresponds to the exchange of more than one 
photon between the electron and proton. If the proton is 
treated as a rigid structure this correction is negligible, the 
difference between an exact phase shift calculation and the 
Born approximation being ^0-25 % [12]. However there 
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exists the possibility of electron-induced polarizability of the 
proton's charge and current cloud [13]. Thus to order e 4 in the 
matrix element the electron may exchange two photons with 
the proton, virtually exciting isobaric states of the proton 
corresponding to resonances in the pion photoproduction 




FIG. II. 3. Feynman diagrams for radiative corrections to the electron 
current in electron-proton scattering. 

cross-section. A diagram contributing to this virtual compton 
scattering by the proton is shown in Fig. II. 4. Although they 
contribute corrections of order a = 1/137 to the main term, 
Fig. 1.2, terms of this type may be significant in the electron 
scattering analysis since in the compton scattering of real 
photons by protons they are responsible for the large resonance 
in the cross-section. Approximate evaluation of this polariza 
bility correction has been carried out using a static 3-3 
resonance model for the proton compton scattering [14] and it 
was established that the form factor analysis according to the 
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Rosenbluth formula is valid to order a = 1/137 up to energies 
ffsf I BeV for electron-proton scattering. However, as the 
laboratory energies for electron-proton scattering studies 
increase into the multi-BeV range the virtual compton 
amplitudes may be expected to play a more significant role. 
At higher energies there will be many important intermediate 
channels in Fig. II. 4 for photon absorption to compete with 
the rapidly decreasing amplitude for the one-photon exchange. 




FIG. II. 4. Diagram for an e 4 contribution corresponding to virtual 
compton scattering. 

This result that compton scattering of the virtual photons 
in the electron-proton scattering plays a minor role at present 
energies whereas the compton scattering of real photons 
exhibits a big resonance [15] can be understood with a simple 
classical picture. Consider a monochromatic light wave, 
a cos cot, incident on a proton which is represented as a charged 
oscillator with resonant frequency co r and damping constant y. 
The constants co r and y correspond to the total energy at 
resonance and the width of the 3-3 resonance in the meson- 
jaucleon interaction. The charge of the proton scatters the 
light wave with no change in phase, the amplitude for charge 
scattering being proportional to cos co. The dipole moment 
induced in the proton by the incident radiation gives rise to 
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dipole radiation which is phase-shifted relative to the incident 
wave and is proportional to 

(co^~o) 2 ) cos cot+yco sin cot 

(^s--ft) a ) a +y a a> a '' ( " ' 

The first term in (2.1) is in phase with and interferes with the 
direct charge-scattering term. It is the term of interest since 
this interference term is, to order a, the polarizability 
correction. However this in-phase term is seen to be odd 
about the resonant frequency co = co r and thus is largely 
cancelled in the integral over the virtual photon spectrum of 
the scattered electron. It is the second or out-of-phase term 
in equation (2.1) which gives rise to the resonant peak in the 
proton compton effect. For electron scattering its contribution 
is enhanced at resonance but nonetheless it is of negligible 
importance since it is of order a 2 relative to the leading or 
direct charge-scattering term. 

Evidently electrons and photons are complementary tools 
for the study of nucleon structure. In compton scattering it is 
the polarizability of the structure that is studied whereas 
electron scattering measures the form factors F-^ and F%. A 
connection between the form factors and polarizability 
emerges only as a detailed consequence of specific meson 
theories and models. Whereas the nucleon is real and on the 
physical mass shell both before and after absorbing the virtual 
photon at the vertex in Fig. II. 1, in the compton process, 
Fig. II. 4, it leaves the mass shell in the intermediate state 
which may contain pions and other strongly coupled states. 

In r6sum6, the above discussion shows that in the range of 
present laboratory measurements, electron-proton scattering 
may be interpreted in terms of nucleon form factors using the 
Eosenbluth formula plus the radiative correction of Schwinger 
to the electron current. Measured form factors from some 
recent experimental data are shown in Fig. II. 5. 

For elastic scattering g 2 <0 and so electron-proton 
scattering measurements yield direct information on the form 
factors for space-like momentum transfers only. In order to 
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measure the form factors for time-like 2 >0, which is of 
interest in itself as well as for dispersion theory discussions of 
the nucleon form factors, as wiU be evident in Chapter III, 
it is necessary to turn to annihilation experiments 

e~+e+^P + P, (2.2) 

in which case j a >[4(Jf)a]. It is hoped that such events will 
soon be within the scope of quantitative laboratory study. At 
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FIG. II.5. Experimental data on proton form factors. For 
4 x 10 26 cm~ 2 what is drawn is -F-(2 2 )/-F(0); what is determined by 
experiment is essentially the combination J?* +2Mlt'$. At q z = 
+4-6 and +9-4 xlO 26 cm- 2 , the ratio b = F^MK has been 
measured by comparing experiments at different energies and angles, 
for fixed g, and found to be 1-2 0-2. The experimental points display 
F&(q 2 )/F%(Q) corresponding to three possible ratios 6 = 1-0, 1-2, and 
1-4. For small g 2 <2 XlO 26 cm" 2 what is drawn is JTffe 2 )/^^); 
the magnetic contribution is negligible in this region. 

this time they are not feasible because of energy and intensity 
limitations. Eeaction (2.2) proceeding left to right requires 
^2 BeV in the centre-of-mass system of the colliding electron 
and positron; proceeding to the left it requires an intense 
anti-proton beam as well as sensitive detectors in order to 
give an observable counting rate since the branching ratio for 
this process is <oc 2 = (1/137) 2 . 
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2 [16]. Since there are no free neutron targets, present in 
formation on neutron form factors is in general less precise 
and complete than that available for the proton. However 
one aspect of the neutron structure, the mean square radius 
of the neutron's charge distribution, (r*)i t N, was the first 
experimental result on nucleon structure. It was originally 
measured in 1947 by scattering of slow neutrons from atoms, 
and it remains today as perhaps the most puzzling of all the 
nucleon structure properties. 

Experimental measurements of <r 2 } ljjv have utilized the 
coherence between the spin and velocity independent part of 
the neutron-electron scattering amplitude with the coherent 
nuclear scattering of the neutron in the atom to which the 
electron is bound. This coherence has been observed in the 
laboratory in three different ways: 

(1) For neutron wavelength comparable with atomic size, 
the interference term exhibits an angular dependence 
described by an atomic form factor since the nuclear 
scattering is isotropic and the superposition of scattered 
waves from the different electrons gives a form factor 
variation. Measurement of an angular asymmetry, 
when combined with known atomic form factors and 
coherent nuclear scattering amplitudes, yields (r*) l)N . 
An experiment of this type was first performed by Fermi 
and Marshall [17] using noble gases so that corrections 
in the analysis of the data from induced polarization 
effects in the closed electronic shells are negligible. 

(2) Similarly this interference term gives rise to an energy, 
or wavelength, variation of the total cross-section, which 
can be measured and from which (r^) l tN can be extracted 
in the same way as above. Such experiments were first 
carried out by Havens, Rabi, and Rainwater [18] in 
liquid lead. 

(3) The coherent scattering amplitude a for forward 
scattering of neutrons in a medium is related to the 
index of refraction of the medium for neutrons by n z = 

7r)2Va, where A is the neutron wavelength and N 
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the number of atoms per unit volume. By measuring 
the critical angle of reflection at the interface between 
two media it is then possible to determine their relative 
index of refraction, and hence coherent scattering 
amplitudes, since cos 6 C njn^ where 6 C is measured 
relative to the interface. Thus for C < and n 2 1 

61 = 

For a precise determination of the coherent neutron- 
electron amplitude two materials are sought for which 
the product density times nuclear scattering amplitude 
is approximately equal but for which the atomic 
numbers are very different. Then the importance of the 
neutron-electron amplitude in determining the critical 
angle is enhanced in (2.3). This method has been 
developed by Hughes and co-workers [19] who used 
bismuth and oxygen as the two substances in forming 
their neutron 'mirrors.' 

Precision measurements by each of these quite different 
techniques, carried out on different materials under different 
physical conditions are consistent with one another. This 
lends great confidence both to the experimental accuracies and 
to their theoretical interpretation in terms of a neutron charge 
radius. In order to see how the neutron charge radius emerges 
from these coherence measurements it is necessary to return 
to the expression given in equation (1.8) for the current 4- 
vector of a neutron and take the non-relativistic limit of its 
interaction with an external electric field representing an 
electron charge cloud. This gives for the neutron's charge 
density 



upon dropping terms of higher order than q 2 and reducing to 
two component Pauli spinors, #, 
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where k N = 1*91 is the neutron moment in units of nuclear 
magnetons. 

The second term in the above expression has been discussed 
and interpreted by Foldy and the reader is referred to his fine 
review [16] of the neutron-electron interaction for more 
details. Present measurements are quoted there with an 
accuracy of 5 % and yield 

,v = [-0126-006] x 10~ 26 cm 2 



or 

(r*) ltN = (0-006) X 10- 20 cm 2 . 

An alternative and more usual expression of these results is 
in terms of a constant potential well of depth F and classical 
electron radius r = cnjm = 2-8 x 10~ 13 cm, which gives rise to 
the same scattering amplitude viz. 



Translating the experimental results into these units, 
F (exptl) = -(4-05 0-2) keV. 

The smallness of the observed charge radius of the neutron, 
<" a >i,tf defies any simple interpretation with meson theory. It 
is to be contrasted with the form factor behaviour of the proton 
given in Fig. II.5. Though the mean square charge radius of 
the proton is not directly measured and may in fact vanish 
also, a smooth extrapolation of the data to # 2 == suggests 
radii <r 2 ) ljp ^ (|~1) x 10~ 26 cm 2 in closer accord with 
intuitive notions based on meson theory. 

The discussion of this section has assumed that the electro 
magnetic interaction between electron and neutron need be 
treated only in Born approximation. Again, as in the previous 
electron-proton scattering analysis, polarizability corrections 
corresponding to the exchange of two photons, can be 
computed and as found earlier, are entirely negligible here. 
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3 [20]. For information on the neutron's moment structure as 
well as for further details of its charge distribution it is 
necessary to consider large momentum transfer experiments. 
To this end an extensive program of high energy electron- 
deuteron elastic and inelastic scattering experiments has been 
pursued at Stanford. In order to interpret these measurements 
in terms of neutron electromagnetic structure, one must have a 




FIG. II. 6. Diagram for inelastic scattering of electrons from deuterons. 
The deuteron is represented with a fat cross-hatched line. 

reliable subtraction procedure to remove all effects of the 
spectator proton in the deuteron. In the absence of a rela- 
tivistic theory of two bound nucleons a precise reduction of 
the observed deuteron cross-sections to a neutron form factor 
is impossible. Uncertainties in the present analyses are dis 
cussed here along with proposals for their eventual elimination. 
Consider first the inelastic scattering process e + D -> 
e f +N + P of Pig. II.6, which Yearian and Hofstadter [20], and 
Sobottka [21] have studied in detail for information on the 
neutron magnetic form factor. This process leads to three 
final particles and thus to a continuous energy spectrum of 
electrons at each scattering angle. Measurements have been 
made of the differential cross-sections for the number 
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<&Q e /, of electrons emerging in a given energy interval 
dE e , and into a solid angle dl e ,, and (by counting all electrons 
emerging in a given solid-angle element) of the integrated 
cross-section, J (d^a D ldE ef dl e ,) dE e , as well. At the large 
momentum transfers required to probe nucleon structure, it is 
a very good approximation to neglect interference between 
scattering amplitudes from the neutron and the proton and to 
express the cross-section as a sum of single-particle cross- 
sections. This is because the interference term is reduced in 
amplitude by the fourier transform of the deuteron wave 
function which is very small at these large momenta. 

There is an approximate sum rule expressing the integrated 
cross-section, i.e. the area under the inelastic curve, as the sum 
of free nucleon cross-sections for the same solid angle [22]: 



This relation is useful for large enough momentum transfers so 
that the neutron contribution to (2.5) is comparable to that of 
the proton [23]. There are, however, significant corrections 
for these conditions which severely limit the reliability with 
which the neutron form factor can be deduced from (2-5) 
[22], [24]. Some of these corrections can be calculated, such 
as kinematical corrections arising from the width of the 
momentum distribution in the deuteron. However, mesonic 
corrections from exchange currents in the deuteron and from 
final state interactions can only be estimated in the absence 
of a satisfactory theory of the deuteron. By relating these 
effects with their important observed role in the photo- 
disintegration of the deuteron at high energies corrections to 
(2.5) of the order of 10% are estimated [25]. 

These same estimates suggest that it is possible to avoid the 
main uncertainties from the meson exchange current correc 
tions by considering the value of the sharply peaked maximum 
in the differential cross-section (see Fig. II. 7). Uncertainties 
which remain in the theoretical interpretation of the peak 
value come from the final state interaction between the 
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outgoing nucleons, to which, the sum rule for the integrated 
cross-section is insensitive, and from the deuteron wave 
function. Durand [22] has extended the original calculations 
of Jankus and given an approximate calculation of the effect 
of final state interaction, using the Signell-Marshak phase 




^ 210 230 250 270 290 310 330 
Electron energy in MeV 

FIG. II.7. Spectrum of inelastically scattered electrons from 
deuterium as a function of the energy of the scattered electron for an 
incident electron energy of 600 MeV and a scattering angle of 135. 

shifts, and found that it serves to reduce the peak value of 
d*a D ldE e , dl e , by ^5% for current laboratory parameters. 
At the same time the peak is broadened so that the area under 
the inelastic curve tends to remain constant. 

The deuteron wave function is, of course, related to the final 
state interaction. Its effect on the peak value, calculated by 
comparing with the matrix element for transitions from the 
deuteron ground state to non-interacting final particles, was 
found to be <5%. It arises as follows: the peak in the 
inelastic electron spectrum at a given angle occurs at an energy 
which is, with neglect of the deuteron binding energy, the same 
energy as for scattering by a free nucleon. Upon integration 
over the undetected recoil nucleons, this peak acquires a major 
contribution from the possibility that one nucleon in the 
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deuteron remains at rest as a spectator, and the other one, 
serving as the 'target' in the matrix element, recoils with the 
full momentum. This particular contribution is proportional 
to the zero-momentum amplitude of the deuteron wave 
function, i.e. to the normalization of the asymptotic wave 



0-9 
0-8 
07 
0-6 

0-5 
0-4 

0-3 
07 

'i 


























5 135 points 
J Small angle points 
"Best fit" for 135 ~ 
ooint 


* 












-i 






















K 
















N 


\ 


















s 


k 


H 


T 















J 


1 


\ 


( 

s; 

i 


\ 


vi 


4 6 


8 


10 12 




14 16 



FIG. II. 8. Experimental measurement of the neutron magnetic 

distribution from inelastic scattering of electrons from deuterium. 

The data is plotted assuming a vanishing neutron charge distribution. 

(Taken from S. Sobottka and R. Hofstadter [21], and see [9].) 

function. In the integration over all nucleon recoils, there are 
in addition contributions to the cross-section which are 
proportional to the fourier transform of the deuteron wave 
function for the various momentum values available to the 
recoiling 'spectator nucleon/ 

The experimental situation [20], [21], [26] at present is such 
that uncertainties in the measurements and in the theoretical 
interpretations are comparable. The latest work of Sobottka 
has been summarized in Figs. II. 8 and II.9. Plotted here are 
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data taken both at large and small angles in order to permit 
in a preliminary way separation of the charge and moment 
structures. Fig. II. 8 is plotted assuming that the neutron 
charge distribution vanishes identically, an assumption con 
sistent with but not required by the vanishing of the mean 
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FIG. II. 9. Ratio of charge to magnetic form factors of the neutron, 
(Taken from S. Sobottka and R. Hofstadter [21], and see [9].) 

square radius (r z )i )N in the neutron-electron experiments 
discussed earlier. Fig. II. 9 places limits on the ratio FI (g 2 ) to 
2MF% (g 2 ) deduced from these same data. 

Note that these data are consistent with a similar magnetic 
moment structure for neutron and proton, as is not surprising 
in view of equation (1.12). Thus in terms of isotopic states 



and J? ! 2(0)/J 7 2(0) =0-06/l-85 ^ 3%. For the largest mea 
sured q 2 , F^(q 2 )^Q'3F^ (0) so that even for a point structure 



28 THEORETICAL APPROXIMATIONS 

for the iso-scalar moment, the corrections should be no more 
than 10 %. This is not the case for the charge form factor, for 
which .Ft(O) = F{(Q) =: J6. 

In order to sharpen the above conclusions it is necessary to 
push the present experiments to greater accuracy or to turn 
to more exotic coincidence arrangements. The limitation in 
the former case comes from the fact that the peak and the area 
of the inelastic spectrum are the sums of two parts, the larger 
of which is due to the electron-proton interaction. The 
neutron form factors emerge with magnified errors when 
expressed as the difference of two experimental numbers, for 
deuteron and for proton scattering, respectively. This is 
evident in the preceding two figures which are computed from 
the electron deuteron inelastic peak which is measured with a 
5 % accuracy and from the electron-proton scattering known 
to better than 10 %. 

Coincidence experiments may be designed to avoid this need 
of subtraction by selecting those nucleon recoils which 
correspond predominantly to neutron scattering, with the 
proton in the deuteron as idle spectator. This possibility can 
be seen easily from the following expression for the transition 
amplitude T* Denote by k the initial electron momentum, 
in the laboratory system, by y (r p T N ), the wave function of 
the deuteron initially at rest, by k r = k q, the final electron 
momentum, by p p and p^, the asymptotic free particle 
momenta for the proton and neutron respectively, and by 
/(r p r A r), a final state nucleon correlation: 



-r 



Here t^i^ r fi ) is the matrix element for the two-body inter 
action of the electron with the ith nucleon and contains all 
spin and normalization factors. Re-expressing T in terms of 
centre-of-mass and relative co-ordinates, with rs=r p T N , 
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and extracting the statement of overall momentum conser 
vation, there results 

T = 
x 

+ [Je-^^ v (x)^x]^ r /*(r)e^-^(r)j. (2.6) 

The two terms above express contributions to the transition 
amplitude arising from electron-proton and electron-neutron 
scattering, respectively. Consider an experiment in which the 
scattered electron is detected with a given final energy* and 
angle corresponding to a large momentum transfer q as is 
desired for the nucleon structure analysis. The first term is 
the dominant contribution for a proton emerging with the full 
momentum transferred from the electron while the neutron 
stays at rest as spectator. This is due to the fact that for 
p_v -> 0, the integral J^ 3 f/ 5is (r)e-' p - v ' r 7/; i7 (r) receives important 
contributions out to r~ radius of the deuteron. On the other 
hand for p# ~Q, the integrand oscillates rapidly for r ^ l/q and 
so it is expected that the ratio 



~ |q.| 2 (radius of deuteron) 2 

The second matrix element dominates for the neutrons 
emerging with the full recoil, and with the proton as the idle 
spectator. Evidently the relative magnitudes of these two 
terms can be controlled by experimental arrangements which 
fix the values of the recoil momenta p p and p A ^ subject to the 
constraint tt+pp+pjv = 0. By fixing |pp|<^|pjv| & * s 
possible to emphasize the neutron-electron scattering ampli 
tude by a large factor relative to that for the proton. One can 
thereby avoid the difficulties of the subtraction procedure as 
applied in the analysis of the peak height of the inelastic 
spectrum. 
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While the coincidence arrangement removes the subtraction 
difficulty by making one of the two terms in (2.6) dominant, 
both terms there exhibit a dependence on the final state 
interactions between the emerging nucleons which is expressed 
approximately (with dependence on exchange and spin 
operators suppressed) by the correlation factor /(r^ r p ). 
Durand [22] has proposed a programme of feasible coincidence 
experiments which to a large extent removes the uncertainty 
due to final state interaction by constructing appropriate 
cross-section ratios. His method yields the ratio of the neutron 
to the proton form factors at a fixed value of 4-momentum 
transfer q* from a measurement of the ratio of the number of 
neutrons to the number of protons emerging in a narrow 
conical region about the direction of the electron 3-momentum 
transfer q. Thus, neutrons emerging in a narrow cone about q[ 
lead to small values of proton recoil, p p = q. p# and are 
predominantly the recoil neutrons from electron-neutron 
scattering. The ratio of the numbers of neutrons to protons 
involves the same matrix element, containing the deuteron 
ground state and final inter-nucleon correlation function, both 
in the numerator and denominator. To the accuracy with 
which one of the two terms in (2.6) dominates the other for the 
experimental parameters of interest, this ratio depends only 
on nucleon form factors and calculable kinematical factors. 
Up to corrections of order 1/-M" 2 in the nucleon interaction 
Hamiltonian, Durand finds for the ratio R of the number of 
neutrons to the number of protons scattered into a given 
conical region about q, corresponding to scattering angle 6 of 
the electron, 



= x [2(Jf (q*) +2M Fl (g 2 )) 2 tan* (J9) +{2 Jf ffi (g 2 )} 2 ] 



X [2(F?(q*)+2MFZ(q*)) z tan 2 ($0) +{2M F$ (g 2 )} 2 ] 
or just the ratio of Rosenbluth cross-sections altered by the 



INELASTIC ELECTROlSr-DEUTERON" SCATTERING 31 

appearance here of 3-momeiita |q| due to the nonrelativistic 
treatment of nucleons with a 2-component Pauli equation. 
Corrections to E arising from interference between the neutron 
and proton terms in (2.6) were calculated by Durand and 
typically found to be very small. For example, at the peak of 
the inelastic spectrum for the scattering of 500 MeV electrons 
through 75, Durand finds, with neglect of final state inter 
action and with an S-wave Hulthen deuteron model that the 
interference term contributes less than f % when integrated 
over a cone of half-angle 30 about the q direction in the 
centre-of-mass system of the outgoing nucleons. When final 
state interactions are included the interference between the 
two amplitudes in (2.6) becomes somewhat more important 
due to the presence of high momentum amplitudes in the 
Fourier transform of f(c N r p ). Durand has shown by 
approximate calculation with the Signell-Marshak phase 
shifts that this correction serves to increase R for the above 
values of experimental parameters by no more than ^5%. 
The effect at larger scattering angles is further reduced [27]. 

Durand's analysis of the ratio R points to a way of obtaining 
greater accuracy in the determination of neutron electro 
magnetic structure from experiments on deuteron targets. 
His programme can be carried through either by comparing 
numbers of protons and neutrons in a cone about q, or by 
comparing numbers of protons in forward relative to backward 
cones in the centre-of-mass system of the emerging nucleons 
about q. The most difficult points in these two experiments 
are the requirements of counting fast neutrons in the first 
arrangement and of counting slow protons in the second. One 
may avoid these by subtracting the cross-section for protons 
in the forward cone about q from the total one d 2 a/dl$ e , dl e , 
for the same q, thereby obtaining directly the cross-section 
integrated over the backward cone to good accuracy. The 
advantage of this subtraction procedure in comparison with 
the one now used for the peak height of the inelastic spectrum 
derives from the facts that it is made within one experiment 
and that the final state interaction f(r) appears in all matrix 
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elements. It is of course absent in the electron-proton 
scattering now used in carrying through a subtraction 
programme. 

In principle it is possible to push the coincidence technique 
one step further and, by an extrapolation procedure, free the 
analysis entirely of corrections due to non-relativistic approxi 
mations, binding, and final state interactions between the two 




FIG. 11.10. A possible decomposition of the Feynman diagram for the 

reaction e+D -*e+N-\-P. The series of parallel lines, labelled fji, 

represents all possible intermediate states coupling the two pa-rts of 

the diagram. 

nucleons. Chew and Low [28] have described a general method 
for such a programme which is based on the possibility of 
analytically continuing the scattering amplitude 



to a pole, the existence of which and the residue at which can 
be plausibly conjectured from field theory. To be specific for 
the present application, consider Fig. II. 6 for the amplitude 
of interest and dissect it as indicated in Kg. 11.10. Here ^ 
represents all possible states coupling the two parts of the 
diagram; for example p represents a neutron, or neutron plus 
pion, etc. (only the strong couplings are included). Consider 
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now the amplitude for this process as a function of the 
invariant square of the four-momentum difference between 
incoming deuteron and emerging 'spectator' proton: 



= (M D ~M P )*-2M D T Pl (2.7) 

where the subscript on the 4-momentum p and on the mass M 
indicates deuteron (D) 9 proton (P), or neutron (N) 9 and where 
T Pi is the proton recoil kinetic energy in the lab system. The 
basic assumption is that the analytic continuation of this 
amplitude as a function of A 2 has a simple pole at A 2 = M% 
for the single particle one-neutron state and a branch cut 
starting at the point A 2 = (M N +m^* corresponding to the 
state of one neutron plus one pion. The practical value of this 
continuation procedure depends on the possibility of extra 
polating the experimental measurements to this pole and, 
hence depends on the proximity of the pole to the physical 
range of scattering parameters and on its separation from other 
singularities in the scattering amplitude. It is of interest then 
to analyse further the kinenaatical questions relating to this 
process and to locate other singularities in the scattering 
amplitude. The general recipe for locating these singularities 
follows from, the work of Mandelstam [29] and may be stated 
simply as follows: Consider first the states which can be 
reached withp^ and p p the incident particles. The outgoing 
proton is replaced by its anti-particle as incident in the 
process, and so 

deuteron+anti-proton -> neutron 

neutron -f pion 
etc. 

There is a pole in the amplitude for A 2 = (p D p P } 2 = M Z N , 
the square of the mass of the one-particle intermediate state 
that can be reached, and branch cuts starting at the squared 
masses of all multi-particle states, the nearest one in this case 
being (M N +m n ) 2 . This gives the spectrum of singularities for 
p D and p p incident. For the remaining singularities consider 
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the crossed reaction with incident deuteron plus anti-neutron, 
(P D ~PN}* = ^- 

deuteron -f anti-neutron -> proton 

proton and pion 
etc. 

There is then a pole for cr 2 = Mp and a branch cut starting at 
(M p+m^ 2 . By relating a 2 and A 2 through the total energy of 
the reaction J57 2 = (p P +p N ) 2 , where E is the total energy of 
the final nucleons in their centre-of-mass system (Pp+p^- = 
0), the singularities in the crossed spectrum are expressed in 
the A 2 plane. It follows from the relation 



that the pole in the crossed spectrum lies at 



and that the branch cut starts at 
A 2 



This spectrum of singularities is readily transcribed in terms 
of cos 6 where 6 is the angle between the direction of 3- 
momentum transfer q. and of the neutron in the centre-of-mass 
system of the emerging nucleons; i.e. cos 6 = q . pV|ft||P&| 
in the system p^ +p$> = 0. (E D denotes the total energy of the 
deuteron in this system.) By overall energy-momentum 
conservation 



s 0). 
The pole at A 2 = M\ corresponds to cos = -fr, i.e. 

^|q. ||p^| 

a EE D -Ml 

nr\o H _ 4%^.^,.. " __ ^^ 1 

cosy o - ^(El-Ml)VW~~-4:M 2 N ) >L 
The branch out starting at (M N +mJ z corresponds to 
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The pole in the crossed spectrum at A 2 = E^+M' 2 N +M 2 D +q 2 
and the corresponding branch cut appear symmetrically at 
cos (TT + O) = cos 6. In the cos 6 plane the analyticity 
domain is as shown in Fig. 11.11. 

The amplitude is analytic in the physical region with an 
isolated pole just beyond the boundary of the real physical 
variables in both the forward and backward cones about q and 



co* 6 plane 



(branch cut) (pole)-1 



(pole) (branch cut) 



1 COS0 COS0J 



FIG. 11.11. Poles and branch cuts in the cos 6 plane for the amplitude 
of the reaction e-j-D -> e+N+P. 

with branch cuts further out. That the pole lies outside the 
physical region is nothing more than a stability statement for 
the deuteron ground state. Equivalently from equation (2.7), 
at the pole A 2 = M^ the recoil kinetic energy of the spectator 
proton in the laboratory system is [30], [29] 

T Plo l(M D -M P -M N ) = -1-1 MeV; 
the corresponding energy for the first branch point lies at 
v lM D ) (Mv-Mp-Mx-mJ ** -70 MeV. 



The singularities in the crossed spectrum are as above with the 
recoil neutron kinetic energy T Ni replacing T Pi . The pole 
corresponding to the one-nucleon intermediate state in Fig. 
11.11 is thus well separated from the other singularities and 
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lies quite close to the edge of the physical region. It is to this 
pole that one extrapolates by detecting protons in the back 
ward cone about q. and studying the cross-section for fixed 4- 
momentum transfer g 2 as the proton kinetic energy T PI goes 
to zero. According to the Chew-Low conjecture [28] the 
residue at this pole is given by the product of the transition 
amplitudes connecting the two groups of particles in Fig. 11.10 
to the intermediate neutron on its mass shell. In this case the 
residue of the pole at A 2 = M 2 N is the product of the amplitude 
for neutron-electron scattering times the amplitude for 
deuteron 4- neutron +proton. Since the amplitude for D < 
N+P can be expressed in terms of the normalization of the 
asymptotic deuteron wave function, the cross-section for Fig. 
11.10 has a second order pole at A 2 = M%- whose residue is a 
known multiple of the electron-neutron scattering cross- 
section. In this way by programming experiments for fixed g 2 , 
the neutron structure can in principle be determined. 

The relation between this extrapolation procedure of Chew 
and Low and the analysis of Durand is now clear. For protons 
emerging in the backward cone about q in the centre-of-mass 
system of the recoiling nucleons. the transition amplitude is 
near the pole at cos 6 = cos , or T PI T PI , corresponding 
to the one-neutron intermediate state in Fig. 11.10 and is 
dominated by the second or neutron scattering term in 
equation (2.6). The closer the experimental parameters 
approach this pole the less important are the contributions 
from the pole in the crossed reaction, corresponding to proton 
scattering, and the less important are the contributions from 
the branch cuts arising from higher mass intermediate states, 
such as neutron plus pion, which correspond to final state 
interactions as expressed by the correlation factor /(r^rp) 
in the form factor integrals in equation (2.6). This result can 
be seen from equation (2.6) as well as from the field theory 
formalism. Note that at the pole, the integral over the 
deuteron ground state wave function diverges since the 
asymptotic exponential decay of y g (r) for r >- oo is com 
pensated by the exponential increase of exp ik p N . r for 
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imaginary k P|JV at the pole T^ = k*/2M = -M MeV. The 
form factor integral then becomes infinite upon integration 
over all space since f(r) -> 1 for r> inter-nucleon force range. 
The final state interaction contributes a finite additive term 
arising from the difference (f(r)-l) which vanishes for large 
r. This term corresponds to the finite contributions from the 
higher branch cuts in Fig. II. 11 and is relatively less and less 
significant as one approaches the pole for the one-nucleon 
state. Hence at the pole, only the asymptotic normalization 
of the deuteron wave function (i.e. the zero momentum 
amplitude in the Fourier transform of the wave function) 
appears in the calculated cross-section and this can in fact 
also be removed by taking a ratio of extrapolations for 
emerging neutrons as well as protons. 

The explicit calculation of the formula for the residue at the 
pole at A 2 = M 2 N in Fig. 11.10 proceeds as follows. The 
matrix element corresponding to this diagram is 

ME = g( Pp )(ME) eN (2.8) 

where p p is the proton recoil momentum in the laboratory 
co-ordinates, and at the pole A 2 = M* N , (ME) eN represents 
just the electron-neutron scattering amplitude on the mass 
shell [31]. The amplitude g(p p ) is the projection of the initial 
deuteron wave function on the recoil proton plane wave state 
and is given by 

(?&(*p> 0, Wfrp* *N, <)) 



Comparing with the earlier discussion of equation (2.6) it is 
seen that equation (2.8) corresponds to the second or neutron 
scattering term and g(p p ) is just the last factor there, expressing 
the amplitude for a spectator proton to emerge with p p from 
the deuteron. Near the pole the internucleon correlation, /(r), 
may be neglected as discussed in the preceding paragraph [32]. 
To a very good approximation g may be expressed in terms 
of the deuteron radius I/a and triplet effective range r by 
, [33], 
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The differential cross-section for scattering of the electron into 
energy-angle interval cQ c , dE e , and of the recoil spectator 
proton into kinetic energy interval dT PI in the laboratory 
system is then, as T PI - %* P , 

_ (da Ne (q*)\ _H N M_p_ 
\k' 



(M N +k,)M P da Ne (q*)\ 

*Q.' J ( } 



where & and Jf Q represent the incident electron energy in the 
laboratory system and rest system of the intermediate target 
neutron, respectively. To complete the kinematics one writes, 
in terms of the square sum of the final electron-plus-neutron 
4-momenta co 2 , 



where, as always, terms proportional to the electron rest mass 
are neglected. The neutron structure is then determined for 
fixed q 2 by the extrapolation 

\ 




The ratio of extrapolations for T Pi -> T P and T Ni > T Nn re 
moves F, the normalization of the asymptotic wave function of 
the deuteron, and fixes the neutron structure in terms of mea 
sured proton form factors. Moreover the absolute magnitude of 
the extrapolation to the proton pole may be studied as a check 
upon the quantitative validity of the entire procedure in the 
polology analysis since there is accurate independent informa 
tion on the electron-proton scattering from the free proton 
target measurements, and F is well known from the low energy 
nucleon-nucleon interactions. 

In r6sum6 of this section, coincidence measurements on the 
inelastic scattering of electrons from deuterons offer the 
possibility of considerably sharpening the determination of 
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neutron electromagnetic structure. The main experimental 
problems they pose are due to the necessity of detecting slow 
protons (or equivalently fast neutrons with high energy 
resolution), and to the small cross-section. For example, the 
ratio in (2.9) of the coincidence cross-section to that for elastic 
electron-neutron scattering into a given solid angle, integrated 
over recoil proton energies of 5 to 20 MeV and over a 10 MeV 
resolution width for final electrons is Hr%> for 650 MeV 
incident electrons. 

In conclusion a variation on the polology theme of Chew and 
Low with an apparent experimental advantage may be 
proposed. It consists of inelastic electron scattering from H 3 
and from He 3 , with the detection, in coincidence, of the 
scattered electron and of the spectator deuteron: 

(a) 

(2.10) 
(b) V } 

The ratio of extrapolations to the pole at a deuteron laboratory 
kinetic energy T Di = -2 MeV for these two reactions again 
measures neutron in terms of proton form factors [34]. These 
processes offer two experimental advantages relative to (2.6) 
and (2.9) for deuteron scattering. (1) In measuring the ratio 
of (a) to (b) in (2.10) the same detection equipment may be 
used since slow deuterons are collected in both cases, and one 
need only change the target material. (2) The counting rate 
may be increased relative to (2.9) since it may not be necessary 
to extrapolate so close to the pole at 2 MeV. This is because 
the competing matrix elements in (2.6) for neutron and proton 
scattering are comparable except near the pole where one of 
them dominates. Here however the competing amplitude 
contributing to the cross-section (2.10) is electron-deuteron 
elastic scattering which, in the energy region of interest for 
nucleon. structure, is considerably reduced owing to the large 
deuteron size. Therefore one is not forced to very low deuteron 
recoil energies to be sure that the process is dominated by 
electron scattering from neutrons or protons respectively in 
the two cases of H 3 and He 3 targets. The first order effects of 
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all but coulomb corrections will disappear in the ratio of 
(2.10a) to (2.10b); the coulomb corrections can be computed 
and will be at most a few per cent due to the small value of 
total charge and to the fact that the relative velocities of all 
particles are large. Once again the applicability of this 
procedure can be directly checked (at least in principle) by 
study of (2.1 Ob) and its analysis in terms of known proton 
structure and of low energy proton deuteron scattering. 

4 [35]. Elastic electron-deuteron scattering may be mentioned 
briefly as a source of corroborative evidence on neutron 
structure. However theoretical attempts to draw quantitative 
conclusions regarding neutron form factors from this process 
again run into the fundamental difficulty that the relativistic 
two-nucleon problem is still unsolved. At present the best one 
can do is to treat the deuteron in the non-relativistic limit as a 
static bound system of two Pauli nucleons and deduce from 
measurements of the elastic cross-section the product of the 
nucleon form factors and the fourier transform of the square 
of the deuteron's wave function. A detailed analysis and 
critique of this approximation has been given by Blanken- 
becler [36]. For not too large scattering angles charge 
scattering dominates and what is measured is 

where 



is commonly called the deuteron ground state form factor. 
The sum of the nucleon charge form factors appears here 
because the two nucleons scatter coherently in an elastic 
process. Note that the combination 



is just the isotopic scalar charge form factor. 

Mclntyre [35] and collaborators have carried through an 
extensive programme of elastic scattering from the deuteron 
and their results are presented in Kg. 11.12. The theoretical 
curves in this figure are drawn with the assumption that the 
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neutron charge distribution vanishes identically for all values 
of g 2 and shows the influence of the proton form factor and of 
various deuteron models in fitting the data. Conversely, were 
bhere a known ^(q 2 ) this experiment would provide F{(q*} 
and hence in conjunction with proton scattering measurements, 
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FIG. 11.12. Elastic eleetron-deuteron scattering data presented as a 
ratio of experimental results to theoretical prediction for a point 
deuteron composed of point nucleons. The upper curve labelled 
proton indicates the proton form factor for comparison. The 
theoretical curves for three different deuteron models are constructed 
with the assumption that the neutron charge distribution vanishes 
and hence represent [# 7 (g r2 )-# T D(<? 2 )] 2 . The argument of the deuteron 
ground state form factor is taken to be the four momentum transfer, 
or the three momentum transfer in the centre of mass system. (Taken 
from J. A. Mclntyre and G. B. Burleson [35]). 



. The following examples indicate the sensitivity of 
bhese results to various assumed F^(q*), consistent with the 
observation that (^ 2 }i^->0. First note that with the 
assumptions that (r*) l>N = 0, (r 2 ) 1)P = (O8/) 2 , and that the 
higher moments of the neutron's charge distribution cancel 
those of the proton, the above curves are decreased to less than 
J- of their drawn values at the highest g-values. Secondly it 
follows from the analysis of Schiff [37] that if a neutron has a 
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charge distribution with (r 2 ) 1)jv = and with no more than 
0-le at a distance of one meson compton wavelength from 
the origin, the cross-section is reduced by more than 60 per 
cent at the high q values indicated in Fig. 11.12. 

It is tempting to proceed further and attempt to draw 
quantitative conclusions concerning the neutron charge 
distribution. That the dangers in so doing are considerable, 
however, is apparent from calculations of Blankenbecler [36] 
and of Bernstein [38]. Their (different) relativistic and ex 
change current corrections to the usual description of the 
deuteron as two statically bound Pauli nucleons produced 
corrections of different signs but of magnitude comparable 
with the difference between the curves drawn in Fig. 11.12 for 
different deuteron models. These corrections limit, at present, 
elastic deuteron scattering as a source of information on 
neutron structure. Once this is known from other measure 
ments such as, perhaps, the inelastic cross-section for electrons 
on deuterons, elastic scattering will provide sensitive measure 
ments of properties of the deuteron system itself. 

Friedman, Kendall and Gram [39] have performed elastic 
scattering experiments designed to provide new information 
on neutron magnetic moment structure. They work at larger 
scattering angles in order to emphasize the magnetic interaction 
and form ratios of cross-sections for different energies and 
angles, corresponding to the same momentum transfer, in 
order to isolate the magnetic from the charge scattering, as 
was discussed earlier for the electron-proton scattering. Fig. 
11.13 summarizes their results. 

These ratios have two main virtues. First of all they are 
sensitive to a difference in the neutron and proton moment 
structures. This is seen from the following elastic scattering 
formula, adapted from Jankus, which treats the deuteron 
non-relativistically : 
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where 

s cos 2 (i0) 1 



Because the anomalous static nucleon moments almost 
completely cancel in the deuteron any difference in their 
structures for # 2 ^ shows up magnified at large angles by its 
sizable contribution to elastic scattering. Also if the experi 
ments are performed at not too high energies so that very 
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FIG. 11.13. Ratio F$(q 2 )[F%(q 2 ) as determined by elastic electron 
deuteron. scattering. The experimental data are analysed here 
according to equation (2.11) with the additional simplifying assump 
tions that F*(q*)/F?(Q) = F?(q*)IF*(Q), which is consistent with 
experiment (see Fig. II.5) and that Fl[(q*) = 0. (We quote here the 
preliminary analysis of Friedman, Kendall, and G-ram^ [39].) 

little recoil energy is transferred to the deuteron and hence the 
three- and four-momentum transfers differ by very little, the 
deuteron wave function and structure disappear in the ratio 
of cross-sections (2.11) at different angles. However there 
still remains the difficulty that the relativistic corrections to 
the description of the deuteron are unknown [25], [36]. They 
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may not be at all insignificant in a process such as this where 
one is working with very small cross-sections. 

In discussing the elastic deuteron scattering, the first Born 
approximation has been used. The validity of this approxi 
mation can be defended here to order 1/137 as in the case of 
the electron-proton scattering. Schiff [40], Lewis [41], and 
Valk and Malenka [42] have studied the dispersion contribution 
associated with the possibility of virtually disintegrating the 
deuteron in the intermediate state in a second Born 
approximation and found this to be small. 

5 [26], [43]. The process of pion production in electron-nucleon 
collisions, 



offers two unique possibilities as a probe of nucleon structure. 
First, as a probe of neutron electromagnetic form factors, it is 
free of complications of the two-nucleon problem which have 
been analysed in previous paragraphs. Secondly, it offers the 
possibility of measuring the electromagnetic structure of the 
pion. In the theoretical analysis developed in the following 
chapter, it becomes apparent that a theoretical discussion of 
nucleon structure hinges upon pion structure; in particular, 
approximate dispersion theoretic calculations yield nucleon 
form factors as integrals involving pion structure. The pion 
form factor is then of great interest, not only in its own right, 
but because theory leads to an approximate relation between 
pion and nucleon form factors which in principle, at least, can 
be directly checked by experiment. 

Before outlining the theoretical campaign from process (2.12) 
to actual form factors, it should be pointed out that with 
present experimental limitations it does not appear to be 
feasible to use electron production of pions to probe the pion 
form factors. However, new data on neutron structure can be 
obtained this way. The first study of electron production of 
pions as a nucleon structure probe was given by Fubini, Nambu 
and Wataghin [44], using the method of dispersion theory. 
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Again to lowest order (1st Born approximation) in the 
electromagnetic coupling, the amplitude for (2.12) can be 
expressed, as seen in Fig. 11.14, as the amplitude for pion 
production by a virtual photon of 4-momentum q^: 



The matrix element of the photon current operator^ is taken 
between initial proton and final nucleon plus pion states and 




FIG. 11.14. Diagram for the process (equation (2.12)). 

differs from the photo meson production matrix element only 
by virtue of the non-vanishing photon g 2 and of the fact that 
its scalar product is taken with the M011er potential of the 
scattered electron, which has longitudinal and scalar as well 
as transverse components. These differences are reflected by 
the presence of electromagnetic form factors in the inhomo- 
geneous terms in the dispersion analysis of the amplitude 
(Psft^lj (Q^pi). The inhomogeneous terms come from the 
one-particle intermediate states shown in Fig. 11.15. These 
terms exhibit poles for the one-particle intermediate states. 
The locations of these poles are at energies for which the inter 
mediate particles have a squared 4-momentum equal to the 
square of their mass; i.e. they are 'on the mass shell/ and the 
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residues at these poles are known coefficients expressing the 
amplitude for pion emission by a nucleon multiplied by the 
electromagnetic form factors. In Fig. II.15(a and b), the pole 
occurs for a 'real' intermediate nucleon and the residue depends 




FIG. 11.15. Diagrams producing poles in the process (equation (2.12)). 
Figure (c) occurs only for charged pion production. 

on the nucleon form factors; in the third diagram the pole 
corresponds to a real charged pion and measures the pion form 
factor [45]. 

Just as in the earlier discussion of inelastic electron- 
deuteron scattering, there exists no rigorous proof of these 
analyticity properties. They follow from the Mandelstam [29] 
conjecture of a two-dimensional spectral representation in 
scattering (in this case of a virtual photon to a pion) and are 
valid at least in perturbation theory. 



ELECTRON PRODUCTION OF PIONS 47 

The work of Fubini, Nambu and Wataghin actually preceded 
that of Mandelstam and was carried out in the framework of the 
dispersion theory method applied to the energy but not the 
momentum transfer in the electron-pion production ampli 
tude. Within this more limited 'one-dimensional 3 dispersion 
theory approach it is easy to see that the first two graphs in 
Fig. 11.15 produce poles at energies (#> 2 +&) 2 = -^ 2 and 
(Pi &) 2 = -^ 2 , respectively; or in an invariant notation, at 
where v = (Pi+Pz) . q/2M and 



v = 

v s = k .q/2M. (2.13) 

Fig. II.ISc gives a pole in the momentum transfer variable 
A 2 = (p2Pi) 2 = (q &) 2 = m* corresponding to exchange of 
the pion on its mass shell. Its contribution to the dispersion 
relations in the energy variable appears in the form of a 
general polynomial in v which the 'one-dimensional' dispersion 
relations themselves leave, in general, arbitrary, but which 
can be tied down to the form predicted by perturbation 
theory. 

No such arbitrariness remains in the Mandelstam repre 
sentation which is a two-dimensional dispersion relation in 
both variables v and A 2 and which specifies the residues at all 
poles in v and A 2 . This has already been seen explicitly in the 
deuteron inelastic scattering discussion which considered both 
the direct and crossed reaction and located the corresponding 
poles and specified their residues. 

The residue in the pole at A 2 = m* 9 which corresponds to 
elastic electron-pion scattering, is evidently proportional to 
F^q 2 ). Since this is unknown, experiments to measure 
nucleon structure must be programmed so as to emphasize the 
role of diagrams (a) and (b) relative to (c). Conversely, in 
principle one can programme to emphasize (c) by extrapolating 
in the Chew-Low manner to its pole at A 2 = m* and thereby 
measuring F w . In either case, the theoretical problem is well 
defined here and free of the two-body complications presented 
by the deuteron. Both in its theoretical and experimental 
aspects, however, there is at present much need of refinement 
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before this can serve as a quantitative tool in electromagnetic 
structure investigations. 

Until now, the experimental investigations of Panofsky 
and Allton [43] have aimed at exploring the nucleon form 
factors under conditions such that the pion scattering term of 




FIG. 11.16. Cross-section for electron production of mesons as a 
function of momentum transfer q 2 and of the relative pion-nucleon 
energy E. The shaded curve at q 2 = is theoretically computed for 
photo -meson production. (The curve and experimental points are 
taken from Panofsky and Allton, [43].) 

Kg. II.15(c) is of secondary importance. These experiments 
detected only the scattered electron energy and angle and the 
runs were programmed by varying momentum transfer q 2 for 
fixed energy of the recoiling pion-nucleon system. 

Fig. 11.16 shows the cross-section behaviour and some 
recent data plotted as a function of two variables, # 2 and E, 
the total pion plus nucleon energy in their centre-of-mass 
system. By riding down the peak of the curve, for fixed E, 
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one measures predominantly the isovector part of the moment 
form factor. This is because ^75% of the cross-section at 
resonance comes from magnetic dipole absorption. The 
matrix element for this is proportional to 



-(<Z)] sin < 

therefore only the form factors vary with q* and the phase 
shift is constant for fixed E. An experimental programme 
which rides down the ledge of Pig. 11.16 near threshold is 
feasible and measures predominantly $-wave ^-production 
(via the Kroll-Ruderman term), and hence 



Theoretical difficulties in this programme centre about the 
problem of obtaining accurate solutions to the dispersion 
relations. Considerable calculational effort is called for and is 
being devoted to surmounting this obstacle [46]. In addition 
an experimental programme has been undertaken to detect 
protons in coincidence with the electrons. This will permit the 
charged and neutral pion processes to be separated experi 
mentally [47] and will also contribute to a sharpening in the 
interpretation possible for this process. 

With the feasibility of coincidence studies comes the 
possibility of an extrapolation procedure to the poles (2.13) of 
Fig. 11.15, which gives the form factors directly as in the 
earlier deuteron discussion. 

For the first diagram the pole corresponding to the real 
intermediate nucleon occurs as in (2.13) at a recoil energy of 
the pion given by 



or 



where t^^ is the pion kinetic energy in the centre-of-mass 
system of the recoil nucleon and pion, (p 2 -fk = 0), and is 
sufficiently far removed from the physical region as to 
discourage application of the extrapolation procedure. 

Extrapolation to the nucleon pole in Pig. II.15(b) is also 
ruled out on practical grounds since it occurs at (p^ Jfc) 2 = M 2 
r t vi = fl^-f m*/2jBf where t vi is the pion kinetic energy in 
the laboratory system [48]. 
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From tills it appears that direct, detailed application of the 
dispersion relations as formulated by Fubini, Nambu, and 
Wataghin is required in order to extract the nucleon form 
factors from electron-pion production data. 

The pole in Fig. II.15(c) lies closer to the physical region 
and invites closer inspection as a possible procedure for 
determining the pion structure. In a manner entirely parallel 
to the discussion given earlier for the deuteron break-up, one 
finds this pole at an angle in the p 2 +k = co-ordinate 
system at 



or at a kinetic energy of nucleon recoil in the laboratory 
system T 2i = m%/2M. Frazer [49] has made a detailed 
statistical analysis of the feasibility of the Chew-Low extra 
polation procedure being applied here to measure F v and has 
found that an order of magnitude improvement beyond the 
present experimental accuracy is required. If achieved, 
direct experimental observation of F v (q 2 ) for spacelike 
momentum transfers will be possible. Actually, as the 
theoretical analysis of Chapter III will show, it is knowledge of 
the pion form factor for time-like # 2 >0 that must be inserted 
into a calculation of the nucleon structure. Such would result 
from the same extrapolation procedure applied to the reaction 



since the virtual photon in this case produces a pair with time- 
like momentum. This experiment unfortunately appears to be 
beyond experimental possibilities at present. 

In conclusion, then, electron-pion production is a desirable 
probe of nucleon structure. It is free of the complications of 
the two-nucleon problem and in addition offers the possibility 
of measuring the pion electromagnetic size which enters into 
dispersion theoretic calculations of the nucleon form factors. 
However at this time both theoretical and experimental 
refinements are necessary. 



CHAPTER III 

CALCULATIONS OF 
NUCLEON FORM FACTORS [50] 

THEOKETICAL attempts to understand the electromagnetic 
structure of nucleons on the basis of meson theory originated, 
naturally enough, with perturbation theory calculations in the 
meson-nucleon coupling strength [51]. The agreement 
achieved was, to say the least, unimpressive, which was not 
surprising in view of the size of the meson-nucleon coupling 
constant. In conjunction with the development of improved 
methods of calculation with meson theory, the next theoretical 
work on the electromagnetic structure was done within the 
framework of the static theory of Chew and Low [52]. Several 
limitations are, of course, inherent to any discussions based on 
a static source theory, so only certain features of the structure 
could be predicted on the basis of the Chew-Low theory. 
These limitations were finally lifted, in principle at least, by 
the advent of the dispersion theory [53]. This language allows 
one to set up the nucleon form factor problem in terms of a set 
of equations of relatively simple structure, from which all the 
earlier calculations can be easily derived as straightforward 
approximations. The language allows, also, a number of 
improvements on the earlier theories, although it is certainly 
true that so far no really unimpeachable solution has been 
found by the dispersion method either. 

1 . The theoretical discussion of the form factors may be begun, 
then, by developing the dispersion theory approach, without 
slighting any of the earlier methods which will appear as 
appropriate limiting cases. 

The dispersion theory is based on certain statements about 
the analyticity properties of the form factors, considered as 
functions of a complex momentum transfer squared. Specifi 
cally, form factors, e.g. F^q 2 ), have been defined for real 
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values of the variable # 2 . One may define functions, which 
will be denoted by F^z), which are the analytic continuations 
into the complex plane of the functions -P 1 (g 2 ) defined on the 
real axis, with 



The properties of the analytic continuation which will be 
needed are 

(i) FI(Z) is analytic in the closed upper half z plane, except 
for a branch cut on the positive real axis, for 2^4m;;. 

(ii) FI(Z) (or in some cases F-^jz) vanishes as z-^co 

anywhere in the upper half plane. 

As yet no rigorous proof of either of these statements exists for 
the case of interest here, though proofs are available for a few 
somewhat similar cases [53], [54]. That the statements are 
true for any order of perturbation theory, however, can be 
proved [55], and it is generally felt that the present lack of a 
complete proof independent of the perturbation expansion 
reflects the difficulty of constructing such a proof, not the 
falseness of statements (i) and (ii). 

In the interest of simplicity the above properties of F^(z) will 
be demonstrated in the lowest order of perturbation theory only. 

To simplify the algebra, assume that protons, mesons and 
photons are all spin zero bosons. The fact that spin is present 
in the actual case changes none of the properties of interest in 
the calculation which will be carried out; it serves merely to 
force the use of two form, factors, and to require a lot of 
spinology. For this case, then, the form factor (remember 
again that there is only one for a spinless proton) is directly 
given by the Feynman matrix element for the reaction 
y+P->P, where the photon is allowed to be virtual. To 
lowest order, this matrix element is calculated from the two 
Feynman diagrams of Fig. III.l. To illustrate the argument 
consider only Fig. III.l (a). Upon writing out the matrix 
element corresponding to this diagram one gets as its 
contribution to the form factor the expression 
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Since it is of interest only to demonstrate the properties of the 
analytic continuation of F into the complex plane, constant 
factors may be ignored; several such factors have been 
dropped here. In accordance with the conventional Feynman 
rules, remember also that the masses in the propagators of 
equation (3.1) are considered to have a small negative 
imaginary part, as a reminder of how to treat the singularities. 




(a) 

Fia. III.l. Feynman diagrams for the lowest order perturbation theory 
contribution of virtual mesons to the nucleon form factors. 

It is worth observing, in addition, that the integral of equation 
(3.1) converges; this is not the case for the actual problem, 
including the correct spins, so that it is there necessary to 
renormalize the vertex to produce a finite answer. Again, 
however, this is a minor complication and of no interest for 
the purposes under discussion. 

The expression (3.1) may now be rewritten by use of the 
familiar identity used in calculating Feynman diagrams: 



00 

= 2 dx-L dx% ( 



For this case, 



so that 
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where 

K === fO ~\~ fCQ 



The last equality follows because^ 2 =M*,p f * = (p+q)* =-M" 2 , 
and therefore 2p . q = g 2 . Performing the integral on d 4 &' 
gives ! 3..^ 





which may be rewritten as 

i i 

F ^~$ xdx j d y(q*+i e )x*y(l-y) l -M*(l-x)*-m& (3 ' 2) 



by use of the substitutions 

y = Xz/x 
x = 1 # 3 . 

The imaginary part from the masses is here exhibited explicitly ; 
since the coefficients of M 2 , m^ and g 2 are all positive, the 
overall effect of the small negative imaginary part in M and 
m^ is to introduce a small positive imaginary part ie into the 
denominator as shown. 

The way to continue F(q 2 ) into the complex plane is now 
clear; one defines a function F(z) by equation (3.2) with 
q z +ie replaced by z. Evidently F(z) is analytic in the entire 
upper half plane, excluding the real axis; evidently F(z)^Q 
as 2->oo; and evidently 

lim F(z) = 



Thus F(z) is the analytic continuation of F(q 2 ) into the 
complex plane, and the region of analyticity is all z for which 
the denominator in equation (3.2) does not vanish. The 
denominator vanishes if 



x*y(l-y) 



CALCULATIONS OF NUCLEON FORM FACTORS 55 

in view of the range of x and y, this can happen for all z which 
are real and greater than 4m^. In conclusion F(z) is analytic 
in the entire plane with the exception of a branch cut on the 
positive real axis from 4ml to - 

This verifies conditions (i) and (ii) in lowest order pertur 
bation theory; as stated already, they may also be verified in 
any finite order. 




FIG. III.2. Contour in the complex plane used in deriving dispersion 
relations for form factors. 

Return now to the general case and accept once and for all 
the properties (i) and (ii). First consider the case in which F(z) 
itself vanishes as z->oo in the upper half plane. 

By Cauchy's theorem, 



w 

fc J CO Z 

c 

where c is the contour shown in Fig. III. 2. Since F(z) vanishes 
at GO, the contribution of the semicircle at infinity is zero. 
Then, since the quantity of interest is 



where the integration is now along the real axis, and the 
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contour has been made to approach it from the upper half 
plane. Using the symbolic identity 

lim - - =p- +i7rd(x) (3.3) 

>o & ^^ % 

where p means principal value, this becomes 



Now taking the real part of this equation gives 





or equivalently 

1 f Im F(q'*} 

W) = ~ 7r 2 *Z' 2 - 
77 J q z q 2 ^e * 

oo 

The limit ->0 is to be understood from now on. 

If there is a region of the real axis in which Im F(q 2 ) = 0, 
the above equation may be used to define the analytic continua 
tion of F(z) into the lower half plane. To be precise, define 



*>-:/ 



ImF(q'*) ,_ 

~-= ff ( ) 



for an arbitrary 2. Now let 2; approach the real axis, from 
above and below, in a region where Im F(q' 2 ) = 0. The same 
value is obtained, so equation (3.4) in fact defines the analytic 
continuation of F(z) into the lower half plane, across this 
region of the real axis. If, on the other hand, z is allowed to 
approach the real axis from above and below at a point where 
Im F does not vanish, one obtains from (3.3) and (3.4), 
1m (F(z) z+i -F(z) z . i ) = 2i Im 



(Here remember that Im F(q 2 ) = Im lim F(z)\. j 

\ lz~^ 2 +* j / 

Thus the analytic function F(z), defined in the entire plane, 
has a branch cut across which the discontinuity of the function 
is 2i Im F(q 2 ), Elsewhere on the real axis F(q 2 ) is real [56]. By 
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assumption this branch cut exists for # 2 >4raf, so that finally, 

*" 



for the physical form factor itself. 

If only F(z)/z is assumed to vanish as 2-^00, the application 
of Cauchy's theorem to this function instead of to F(z) itself, 
together with a trivial modification of the previous discussion 
to allow for the pole introduced at z = 0, produces a slightly 
different equation, namely: 



Both of these equations are called dispersion relations; the 
first is said to be an unsubtracted, the second, a once 
subtracted dispersion relation. 

2. The dispersion relations are obviously not sufficient in 
themselves to be used as the basis for any calculations. It is 
necessary to find in addition an expression for Im F(g*), or what 
is equivalent, for the discontinuity of F(z) across the branch 
cut. It is also necessary to decide which kind of dispersion 
relation, subtracted or unsubtracted, to use for any given F, 
and if a subtracted one is used, to decide on the value of -F(0). 
To obtain the imaginary parts and to discuss the question of 
subtractions it is convenient to return to the original definition 
of the form factors. They were given by (neglecting isotopic 
spin for the moment) 



and may be rewritten in a different form [57]: 




Here the + symbol denotes the time ordering of the two field 
operators. The differential operator may be moved through 
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the time ordering with the effect of adding an extra term. 
coming from the operation of the time derivative on the step 
functions implied by the time ordering. Explicitly 



Q 
where n(x Q ) = ( Q ^ <Q 

For the standard meson theory, the commutator may be 
evaluated directly [58]; 



where e is the unrenormalized electric charge; e Q Z^ = e, the 
physical renormalized charge. Thus 



^W^- (3 ' 7) 

Several alternatives suggest themselves for the forms of F^ 
and F 2 . Perturbation theory [51] indicates that F^q 2 ) 
diverges logarithmically as g 2 -*co while F z (q 2 )-+Q; if this 
were the case, a subtracted dispersion relation would be 
required for F^ while an unsubtracted one would do for jP 2 . 
As a result of the term proportional to e Q in equation (3.7), 
which contributes an additive constant e^Z\ to F v another 
possibility is that the integral in equation (3.7) vanishes as 
2 2 -^co 5 so that jFj-^/Zj. In this event a subtracted 
relation is again needed for F lt but e^Z\ could be computed 
from it. In either event, the use of a subtracted relation for F^ 
introduces no difficulty since the necessary constant ^(0) is 
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just the electron charge e. It could also be the case that 
j?\->0, allowing the use of the unsubtracted form (3.5) for F. 
In this case one might hope to be able to get some information 
about the charge, which now appears as an integral over Im F^ 
[59]. Finally, there is no proof that either F : or F% or both do 
not diverge so badly as # 2 ->co that dispersion relations with 
even more than one subtraction must be used. In any case, it 
will be assumed at this point that F 2 satisfies an unsubtracted, 
and FI a subtracted, dispersion relation; a subtracted 
relation for F v it should be noted, does not preclude the 
stronger statement that ^ I 1 -^0 with the consequent validity 
of an unsubtracted relation. 

Aside from the 6 term in equation (3.7), then, evaluation of 
the integral results in 



(38) 



In deriving this a complete set of Heisenberg states of energy 
E n and momentum P n , denoted \n), has been inserted. The 
choice of boundary conditions on these states is a matter of 
convenience. From equation (3.8) an explicit formula for the 
imaginary parts of the form factors can now be obtained. 
First consider the second term. The denominator can vanish 
only if 

PI = ^ = M *. 

It is easy to see that this condition can be satisfied only if the 
state n is a one-nucleon state; if it is such a state, however, 
(n\j N (0)\Q) = 0, so the second term is never singular. To see 
this, note that if (n\ is one nucleon of momentum n, then 



= 0, 
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by the Dirac equation. It therefore contributes nothing to the 
discontinuity of the form factors. 

In the first term, a singularity can occur if 

PI = (p 1 '-p) 2 = i\ 

i.e. only for positive real values of q 2 . The discontinuity in this 
term, then, occurs only if g 2 ^4m^ (assuming the state n = two 
pions to be the state of lightest mass which appears), which 
coincides with the position of the branch cut obtained from the 
perturbation diagram in which the pion absorbs the photon. 
The value of the discontinuity is 



x 



Since the discontinuity in F equals 2i Im F, there results 
Im F+i Im F^a^u,) 

) . (3.9) 



The discontinuities, or imaginary parts, then, receive contri 
butions from various intermediate states n. It is convenient to 
symbolize these contributions pictorially as shown in Fig. 
III. 3. Evidently the delta function in (3.9) implies that a 
given state contributes only for q 2 greater than the rest mass 
of that state. If it is desired, the form factors may, as men 
tioned in Chapter I, be defined in terms of the production 
process y~+N+N; this allows the deduction again of 
equation (3.9) with u v replaced by the antiparticle spinor v-. 
The corresponding diagrams are then altered by having the 
incident nucleon line turned around, to become an outgoing 
antinucleon of momentum p, so that q 2 = (p'+p) 2 . This 
description will generally be used in what follows. 

A final complication is to include both neutrons and 
protons by introducing the isotopic spin and isotopic scalar 
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and vector form factors. How this is done has already been 
discussed in the first chapter; therefore, it -will suffice here 
to say that it is only necessary to replace Im _F 1)2 by 
Im -? 7 2,i+ T 3 Im F in equation (3.9) if charge independence 
is assumed. 




i _____ 



FIG. III. 3. Diagrams representing the contributions of various 

intermediate states to the imaginary parts of the nucleon form 

factors as given in equation (3.9). (a) is a general state n; (b) is the 

two-pion, (c) the three-pion, and (d) the nucleon-antinucleon state. 

The problem is set up, then, by combining equation (3.9) 
with the four dispersion relations 



-t> - i e+ ! f Im ^g' 2 ) 

(2 ) - ?e+ v j q >2 {q ' 2 _ q 2_ ie 



(3.10) 



One additional observation should be made about the above 



equations, and that is, that Im 



and Im F s z vanish for 
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thus the thresholds on the dispersion integrals for 
the scalar form factors should be changed from 4m^ to 9m^. 
To see this, recall that the contribution to Im F from any 
given intermediate state is given by equation (3.9). Consider 
the lightest intermediate states, which are those containing 
only 77 mesons, in that equation. The contribution of a state 
of nir mesons is proportional to the amplitude that a photon 
produces n pions, times the amplitude that n pions are 
absorbed by a nucleon line. It will be proved that this state 
contributes only to the isotopic vector from factors if n is even, 
and only to the isotopic scalar form factors if n is odd. 

Since the state of n mesons is produced by a photon, it must 
have zero total charge and be odd under charge conjugation. 
The pions may be described as created by the Hermitian 
fields q> l9 <p 2 9?3 defined by 

<Pi ^ 



1 
/2 



in terms of the usual charged and neutral fields <p and cp Q [60]. 
Now <p creates negative charge, 9?* creates positive charge, and 
<59 zero charge. Then both (p l and <p 2 create a linear combina 
tion of positive and negative charge. In order that the total 
charge of the state of n mesons be zero, therefore, the number 
of particles of type 1 (created by 9^) plus the number of type 2 
must be even. Under charge conjugation, that is, the inter 
change of positive and negative charge, <p-*<p* and <PQ-><P Q , 
and so (p^-^cp^ q> 2 -^ <p 2 , 9?s->9V If the state is to be odd 
under charge conjugation, then, the number of 2's must be odd, 
and therefore so must the number of 1's. Finally, mesons of 
type 1, 2 and 3 absorb on the nucleon line with isotopic 
operators T 15 r 2 and r 3 , respectively. Therefore the total 
isotopic operator for the ^-meson state is r^T^rl'** = 
T n-odd _ r ^ n j g evenj an( j = i if n i s odd. Thus n even or odd 
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leads to an isotopic vector or scalar form factor, and con 
sequently the lightest state which contributes to a scalar form 
factor is the three-pion state, with a threshold 9mf . 

The dispersion relations may be used to derive expressions 
for, e.g. the 'radii' of the charge and current distributions. 
These, it will be remembered, were defined by 

z =- (3 - n) 



If F satisfies a dispersion relation, either subtracted or 
unsubtracted, 



If in addition F were to satisfy an unsubtracted dispersion 
relation, one could write 



Now, if Im F(q' 2 ) is everywhere positive, equations (3.11, 3.12, 
and 3.13) provide a bound on the radius, given by 



f r the vector radii 



radii. 

Since there is no theoretical reason whatever to believe that 
Im F(q'*) does not change sign, however, these bounds are 
useless, and the radii will not be discussed further. However, 
if Im F(q' 2 ) were positive definite, one could obtain predictions 
about the behaviour of F itself which may be tested by 
experiment. For example, if one computes the curvature of 
l/F [61], one obtains 

d* 1 /I 



( { ImF Y_ f lmF f ImJ 
U (g /2 ~g 2 ) 2 / Jq' 2 -q*] (q' 2 -q* 



Im 
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If Im jP(g' 2 )>0 everywhere, and if g 2 <0, which is the region 
in which F(q*) is measured, the denominator here is positive. 
The numerator may be rewritten as 



Im 



(_L___L_ 



which, under the same assumptions, is <0. Thus, a positive 
definite imaginary part for a form factor satisfying an unsub- 
tracted dispersion relation requires the reciprocal of the form 
factor to be concave downward in the region of spacelike q 2 . 
Present experimental evidence is still too scanty to test the 
validity of this prediction. A somewhat easier comparison 
with experiment seems to be afforded by the curvature of 
log F(q 2 ) rather than of l/F(q 2 ) [62]; namely 



/ 1 \ " / F'\ 2 

(tag *)'--*(*) +( T ), 



which is positive under the above assumptions. Additional 
experimental information on the charge and moment form 
factors separately is required to test either of these results. 

3. Up to this point no approximations have been necessary 
(provided one accepts the assumed analyticity properties of 
the form factors) ; equation (3.9), together with the dispersion 
relations, are a formally correct set of equations from which 
the form factors can, in principle, be determined. A practical 
calculation of the form factors, however, will of course require 
some approximations in the use of equation (3.9) to determine 
the imaginary parts. It is convenient to begin the second i.e. 
the approximate phase of the theoretical discussions by 
indicating how the perturbation theory calculations follow 
from the general formalism. 

Perturbation theory takes the view that only terms with the 
smallest number of powers of g, the meson nucleon coupling 
constant, will be of importance. It requires only a moment's 
thought to convince oneself that the only intermediate states 
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in equation (3.9) which contribute to order g 2 are the two-pion 
and the nucleon-antinucleon states. All other states are at 
least of order g*. Pig. III. 4 illustrates this. To evaluate the 
Im F's to order g 2 , then, it is necessary to drop all but these 
two intermediate states in equation (3.9) and to replace the 








FIG. III. 4. (a) and (b) show the two Feynman graphs contributing to 

the form factors to order g 2 , (c) and (d) are the corresponding 

diagrams indicating which intermediate states the same graphs 

contribute to in the dispersion relation language. 

matrix elements which appear in the states retained by their 
lowest-order values. Thus, one should use, e.g. for the 
contribution of the two-pion intermediate state, the expression 
i Im 



(3.u) 
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(where g l3 g 2 and a l9 cr 2 are the momenta and charges of the two 
intermediate TT'S), and write, to order <7 2 , 






A small amount of algebra, together with the identity 



, 



permits a more convenient form for equation (3.14): 
Ov|r 3 Im -Fjfov+irs Im F^ v q v v^) 27T 



ieg 2 



(3.15) 

It is interesting to observe that this is nothing more than the 
standard Feynman amplitude for the diagram of Fig. (III.4(a)), 
multiplied by 2-n- 2 and with the Feynman propagators of the 
intermediate mesons replaced according to 



A similar result appears if one calculates to order g* the 
nucleon-antinucleon state contribution in equation (3.9). 

It is a definitely non-trivial task to show directly that this 
modification in the Feynman matrix element is what is 
required to give its imaginary part. For a Feynman matrix 
element with a single propagator, the analogous identification 
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is simple, because the small negative imaginary part in the 
mass in each propagator produces the identity 






from which one sees that the imaginary part is clearly obtained 
by the replacement 






However no such simple argument works in the present 
example, where three propagators appear. 

If now these expressions are used to evaluate Im F^l to 
order # 2 , and the resulting functions are inserted into the 
dispersion relations, the usual perturbation theory is obtained 
[51]. It should, perhaps, be remarked that it is the usual 
renormalized perturbation theory which is obtained; this 
approach through the dispersion relations is automatically 
finite at every stage and expresses the answer directly in 
terms of renormalized quantities. 

It may be worthwhile to indicate how this derivation of the 
perturbation result goes in a bit more detail. In equation 
(3.15), use the fact that 



and then write 

K'|r5( 

Next observe that 

2 T ffi T Oo(^ ll 5 ffa 2~^ 1 2<'<r a l) = ^Tg T^ = 

CT l r 2 

Inserting these relations in equation (3.15) yields 
Im Fl+i Im 






^^ 
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where 

Q = Itei-tfa) 



P = 

Before evaluating the integral over d*Q it is convenient to 
carry out the spinology, and identify Im J^ and Im F 2 
separately. This is easily done and results in 



Im J% 2 ) 

i r(p.p) 2 -P 2 Q 2 , , 

X ( p_Q)2_lf2L 2P2 + 

Im P( ff 2 ) 



x 



Finally, evaluating these integrals gives at last 



and 



(tan- 1 2/ - (l - -i- tan- 1 y ) } (3.16) 
I tfoV Z/o // 

^ _2g, r g JV r 
and 



x 
where 



These expressions are rather complicated, and the derivation 
of them has been presented in some detail because the 
techniques of calculation are typical of the type of algebra 
necessary in using these dispersion relations, not because these 
particular results are of any great significance. 
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As indicated above, the results quoted here are only for the 
meson current contribution to the perturbation result in Fig. 
III.4(c). Similar calculations will produce the answer for the 
nucleon current part, Fig. III.4(d). 

This perturbation calculation gives the following values for, 
e.g., the static anomalous magnetic moments [51]. 



The agreement with the experimental values is not startling, 
as is only to be expected from a solution obtained by an 
expansion in powers of 15. The difficulty comes from a large 
isotopic scalar contribution by the nucleon current; i.e. from 
the second of the two perturbation diagrams, Fig. III.4(b). 
Note that the meson current contributes only to the isotopic 
vector part, its second-order contribution being l-6(e/2Jf) for 
2 /477 = 15. This result is increased by less than 10% when 
fourth order terms in the meson-nucleon coupling constant are 
included. On the other hand, the fourth order contribution 
from the nucleon current greatly alters the second order 
results of l-55(e/2Jf) and +-57(e/2Jf) to and 4-4(e/2Jf) 
for the isoscalar and isovector parts, respectively. It is there 
fore necessary to look for an improved approximation method. 

4. The dispersion theory approach, as outlined above, is not a 
foundation which suggests that a perturbation expansion is a 
very natural kind of approximation to try. A more plausible 
scheme might be as follows. Observe that a particular inter 
mediate state in equation (3.9) contributes only with a threshold 
equal to the square of the total rest mass of that state. Thus, 
a heavy state appears only for large q' 2 in the Im jP's; in fact, 
for 4ra 7r 2<#' 2 <16 m%, the Im F's are given exactly by only the 
2?r intermediate state. The heavy state contributions in the 
dispersion integrals may therefore be small, due to the 
denominator (g' 2 q 2 ) in the dispersion integral, particularly 
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if one is only interested in evaluating the form factors at the 
negative values of g 2 for which they are measured. The 
dispersion relation approach is based on the hope, then, that 
it is a sensible approximation to drop the effects of heavy 
intermediate states in equation (3.9). It may be mentioned 
that the subtracted dispersion relation has an extra power of 
g' 2 in the denominator which provides an additional weighting 
against heavy state contributions to the imaginary parts. 
One would therefore expect to have more confidence in this 
approximation scheme if used in conjunction with the 
subtracted form of the dispersion relation. 

An argument such as the above, is, of course, hardly 
convincing; nevertheless it seems more reasonable than the 
perturbation theory. Thus the next try might be to say that 
only the lightest state counts; that is the 2-jr state for the 
vector form factors, and the STT state for the scalar form 
factors. Consider first the isotopic vector case; so return to 
equation (3.9), drop all but the 2rr state, but now do not 
replace the remaining amplitudes by perturbation theory. 
When this is done, and the limit M >co (i.e. the static 
nucleon limit) is taken, the results of calculations made with 
the Chew-Low form of the static theory are reproduced [52], 
[53]. Since the two-pion state will be discussed shortly, there 
is no point in going into the algebraic details here. Suffice it to 
say that the values of the static anomalous moments obtained 
from this theory are 



(Since the two-pion state is a pure isotopic vector, the proton 
and neutron moments are of course equal and opposite; the 
static theory is totally incapable of discussing the isoscalar 
moments.) These numbers result for a static theory coupling 
constant 

/* = o-08 = 

and a cutoff 
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This result is somewhat more encouraging than that of 
perturbation theory; however, talking in terms of the static- 
theory is always rather unsatisfactory since these results all 
depend sensitively on cutoffs and are not gauge invariant, 
among other shortcomings; in addition in this application the 
static theory has been shown to introduce serious kinematical 
errors [63]. 

The previous discussion has sketched briefly how the results 
of both perturbation theory and the static meson theory may 
be derived from the general formalism of dispersion relations. 
The natural question to be discussed next, then, is how to 
improve on these approximations. 

The partial success of the static theory, within its obvious 
limitations, supports the hope that the approximation of 
sticking to the lightest intermediate state, even in the rela- 
tivistic case, may not be too absurd, at least if one confines 
oneself to properties of the form factors for small momentum 
transfers. This hope may be reinforced by looking in some 
detail at a particular heavy intermediate state, that of a 
nucleon-antinucleon pair, and showing that it is unlikely that 
this state, at least, contributes in a very important way in 
determining the low momentum transfer form factors. The 
reason that limitations may be placed on the effect of the 
nucleon-antinucleon state is more easily illustrated by 
discussing the case where the spin of the particles is neglected 
The essential features of the argument are all present in this 
case; the presence of spin only serves to complicate the 
algebra. 

It is obvious, from equation (3.9) and the interpretation of 
this equation in terms of pictures, such as shown in Fig. III. 3, 
that the contribution to the imaginary parts of the form 
factors from the NJ$ state is proportional to the nucleon form 
factor itself times the NN scattering amplitude. Since the 
threshold of this contribution is #' 2 = 4=M 2 , the total centre-of- 
mass energy of the NN scattering process is always greater 
than or equal to 2M, so the amplitude appearing is actually 
a physically observable quantity [64]. In addition, the initial 
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nucleon and antinucleon of this scattering have been produced 
by a photon, which has spin one. So, in their centre-of-mass 
system, only the amplitude for scattering in a state odd under 
charge conjugation and of total angular momentum one will 
be relevant. For the case of the actual spins, this limits the 
scattering to be in the SS l and SI^ states; for the illustration 
problem, where the nucleons are given spin zero, the scattering 
is then in the P state. Thus one may write, in the fictitious 
case, 

Im F(q' 2 ) N $ = F(q'*)* sin de#, (3.16) 

where <5 is the P wave NN scattering phase shift, evaluated 
as a total centre-of-mass energy given by V(#' 2 )- Two 
comments are necessary concerning this equation: first, only 
one form factor is involved. This is because spin zero nucleons 
have only a charge form factor and no moment form factor. 
Second, the right-hand side involves F*, not F. This may be 
understood by referring back to equation (3.9), and observing 
that if the matrix element of the nucleon current is to represent 
NN scattering, the intermediate nucleon-antinucleon state in 
equation 3.9 must be a state with outgoing wave boundary 
conditions [65]. But this means that the photon current 
matrix element has as its final state an outgoing NN pair 
instead of an ingoing one ; since the form factors were originally 
defined in terms of the matrix element of the photon current 
with ingoing boundary conditions, F* rather than F appears 
here. To be specific, the relation of the matrix element 



where the final state has ingoing boundary conditions, to the 
matrix element 



in which the final state differs only in that it has outgoing 
boundary conditions, may be obtained from time reversal 
invariance. The form factor for spin zero nucleons is defined 
in terms of the first matrix element by 



(p'-p) f F(q*). 
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Then, to evaluate the second matrix element, note 



where J^(0) r is the time reverse Heisenberg electromagnetic 
current [66], and 

tr\\ J*A(0) if P * s the time component 

3 [A\ ) T 

~3v(ty if j is a space component. 
Therefore, 

0}* 



where 

, , , __-^ __ (p f p)n if ^ is the time component 

~(p' $>){* if /^ is a space component. 
Thus finally 



that is, the only effect of using outgoing states instead of 
ingoing states is to replace F by _F*. 

Consider the contribution, now, of this NN state to some 
low-energy property of the form factor, e.g. the radius. This is 
given by 



6 1 

4M 2 



Now Im F N $(q' 2 ) is given by equation (3.16); hence 



If F(q' 2 ) is assumed bounded, e.g. by its value at zero, there 
results 



Such a limitation is called a 'unitarity limitation/ because 
it has its origin in the statement that |sin <5e z<5 |<L In the 



74 CALCULATIONS OF 3STXJCLEON FORM FACTORS 

actual case of spin J nucleons, similar limitations may be 
derived effectively by the .method indicated, though the 
derivation is complicated by a certain amount of algebra. 
These limitations exist for the radii, static moments, etc. If 
they are valid, they imply, owing to the large mass of the 
nucleon, that the NN state contributes little to the low-energy 
properties of the form factor [67]. 

Several things should be emphasized, however, to avoid 
undue joy at this result. First, limitations can be imposed 
only on the moments of F, so nothing is said about the effect 
of the N$ state on the behaviour of F at higher momentum 
transfers, say comparable to M 2 . Second, to get the limitations, 
it was necessary to assume that | F\ itself was bounded; in fact 
bounded by a value comparable to its value at the origin. This 
assumption is not necessarily correct. For example, if one 
rewrites equation (3.16) in the form (taking <5 to be real, for 
the moment) 

Im F = tan <5 Re F, 

it becomes less obvious that there should be a bound. Only if 
we say Re F becomes zero when <5 goes through |TT do we get 
back the bound on Im F. But it is not necessarily more 
reasonable to say |JP| is bounded and therefore Re F->0 than 
that Re F is finite at resonance, which implies j F\ is unbounded. 
All in all, then, these unitarity limitations must be taken with 
at least a small grain of salt [68]. 

In any case, for one reason or another, one may try to adopt 
the optimistic view that only the two-pion intermediate state 
is important. This, of course, restricts the discussion to the iso- 
topic vector form factors ; the scalar ones will be discussed later. 

5. Analytically, the two-pion contribution to the imaginary 
parts of the form factors is completely given by equation (3.14), 
where the two matrix elements 



and 
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are treated exactly and not replaced by perturbation theory. 
The corresponding diagram representation Fig. III.3(b) gives 
a pictorial idea of what is involved in calculating this contri 
bution. As the picture indicates, the desired contribution will 
involve the amplitudes for the processes 

V 



The first of these is the electromagnetic form factor of the 
pion; to evaluate this one could pursue the same type of 
attack that is being attempted on the nucleon form factors. 
Define, to be specific, the pion form factor by 



. (3.18) 

This is the same form that was used in the simplified discussion 
of the unitarity limitations imposed on the nucleon-anti- 
nucleon state, when the nucleons were taken to be spinless. It 
is perhaps worthwhile to indicate how invariance arguments 
show this form to be the correct one. 

The matrix element {q^a^ 22 2~ > |4( )|0> is a Lorentz 
vector, and depends only on the 4-vectors q^ and q 2 . It is 
therefore a linear combination of (q^+q 2 ) 9 and ({fr g 2 )- If 
YJ/Z = is imposed, the combination (jh+g^) is clearly 
impossible. From the two vectors one can form three scalars, 
and with the two mass conditions q = gf = m* only one 
independent scalar remains, which may be taken to be 
(21+22) 2 - -^ s f ar as th-e isotopic dependence is concerned, it has 
already been shown that the two pions must be in a T 1 
state; hence the factor 



The remaining factors of i/-\/(4a) ]L co 2 ) are put in explicitly so 
that F n (0) = e, the physical meson charge. 

One may then write a dispersion relation for F v , and derive 
an expression analogous to equation (3.9) for Im F v9 which one 
may attempt to evaluate in some approximation. Now the 
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approximation which is being used for the micleon form factor 
includes only the intermediate state of the lowest mass. It 
therefore seems reasonable to try the same approximation for 
the pion form factor. This intermediate state is just the two- 
pion state. Thus, in this approximation, one has exactly the 
situation which was used as an example to illustrate the use of 
the unitarity conditions in limiting the nucleon-antinucleon 
contributions to the imaginary part of the nucleon form 
factors; it is therefore possible to write immediately, to the 
accuracy of keeping just the two-pion state, 

Im ^(g /2 ) = Fl(q'*) sin de id (3.19) 

where <5 is the P-wave TTTT scattering phase shift. 

One comment is worth making here before proceeding. 
Equation (3.19), together with the fact that Im F v is of course a 
real number, implies that F^ has the phase <5. Now equation 
(3.19) is just an approximation, so in general this statement is 
not exact; however (3.19) is correct for g' 2 ^16m^ (below 
threshold for the next possible state, that of four pions), so in 
the region 4raf<g' 2 <16m^ it is in fact a rigorous result that 
the phase of F^is d. 

Coupling (3.19) with the dispersion relation for the pion 
form factor gives the equation 



q 2 
-+ |r J 



2 tan (5 Re F'* 



This equation may in fact be solved exactly for F v in terms of 
6 [69]. In order to minimize the algebra while explicitly 
demonstrating the method of solution it is convenient to work 
with the unsubtracted form of the dispersion relation. Then 

the equation is ^ 

i r+.a-r, &{* r \ T?^ F ( x '\ 

dx' 



with the restriction ^(0) = e. Define 

? tan 6(x') Re F v (x') 



X Z 



dx'i 
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thus FJ(z] is the analytic continuation of FJ(x) into the 
complex plane. Clearly if 

F v (x) = lim F,(z) 

z-*xie 

then 

F,(x) = F,(x)+ 9 
also 

Re F v (x) + = Re F v (x)_ = Re F,(x). 

Now, assuming again that <5 is real, one can write as before 

Im F w (x) = tan d(x) Re F v (x). 
Thus 

^(^i - (lt tan d(x)) Re ^(a?), 
so that 

F v (x) + __ I +i tan d 



__ 

X()- "" 1 -i tan (5 * 
Now define 

L(z) = log F v (z). 

Then L(z) is analytic, since FJ(z] is (except that L(z) has a 
branch cut on the positive real axis from X Q to oo), provided 
F n (z) has no zeros. The discontinuity across this branch cut is 

L + (x)-L_(x) = log = Xd(x). 



Then, using a subtracted dispersion relation for L(z), assuming 
L(z)/z->Q as 



or, finally, replacing x by q z again, 



snce 



This same result woidd not be much more difficult to obtain 
even if the subtracted form of the dispersion relation for F^ 
were retained. F v could consequently be considered known, 
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if the P-wave TTTT phase shift were known at all energies. What 
can be said about this phase shift will be discussed later. 

6. To complete the calculation of the two-pion contribution to 
the Im .F^g s, the annihilation amplitude for the reaction 



is required. This amplitude is, analytically, the matrix 
element 



By the usual invariance arguments, this may be expressed in 
terms of four unknown scalar functions of the two independent 
scalar variables which can be formed from the three vectors 
p', g^andtfg, viz: 



, .-*!> 
^ + 






Here A and B are functions of, for example, # 2 = 
fei+^2) 2 = 4 -^ 2 and (P f ~ ^Y= W*. The use of this ex 
pression, together with the definition of the pion form factor 
from equation (3.18) allows the two-pion contribution to the 
Im F's to be given explicitly in terms of A , B and F n by 
equation (3.14). The algebra necessary to accomplish this is 
identical with that used in the derivation of the perturbation 
theory from equation (3.14) and will not be repeated here. 
The result is [53] 

Eta J?J(a) = 
Im 
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where 



P.QIP*A~(W*,q*), 



((3(P . Q)*- 
Here P, Q, / are defined as before, namely: 



so TF 2 = (P+G) 2 . 

One approach to the evaluation of J^g 2 ) and J 2 (g 2 ) which 
has been attempted is based on the following observation. 
The annihilation amplitude in question, in terms of which 
J. , B and hence J 1 and J 2 are defined, may be re-expressed 
as 



The right hand side of this equation is nothing more than the 
matrix element for pion-nucleon scattering, from an initial 
state consisting of a nucleon ~p and a meson g 1 cr 1 to a final 
state of a nucleon p' and a meson ^ 2 or 2 . Since the final 
meson, for example, has negative energy, this is not a physi 
cally observable scattering process. Nevertheless, it is in 
principle possible, using dispersion relations for meson 
nucleon scattering, to extrapolate the four functions J. , B^ 
from their measured values in the physical region for meson- 
nucleon scattering to the range of variables required in the 
present calculation. This extrapolation has not, however, 
yielded very good results in practice, as there are a number of 
ambiguities and divergences connected with the method [53]. 
This approach will therefore not be pursued further here. 
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More fruitful results can be obtained on the basis of making 
another assumption, which is that the functions J^ and J 2 
satisfy dispersion relations in g 2 . This assumption is the same 
as the statement that each partial wave amplitude in a 
scattering process satisfies a dispersion relation. It is im 
possible to justify this assumption rigorously, as it was 
impossible to justify the dispersion relations for the form 
factors themselves. However, the assumption may be made 
quite plausible through the Mandelstam representation [70]. 
The details of the use of Mandelstam's representation to 
justify, at least in part, the necessary properties of J^ and J 2 is 
too complicated to be entered into here since there is no 
rigorous proof anyway one might as well just accept the 
assumption in the form most convenient for present purposes 
[71]. 

Therefore, the next move will be to investigate the con 
sequences of the statement that J l and J 2 satisfy dispersion 
relations; namely ^ 

1 



and similarly for J 2 . ~ OT 

To proceed from this point it is easiest to return to the 
interpretation of J x and J 2 in terms of the annihilation process 
7T+-77 ^ZV+JV"; thus they may be written pictorially as in 
Tig. III.5(a). 

The contribution of the two-pion intermediate state, shown 
in Fig. III.5(b) is as before 

Im J*! = J x sin 6e id 

Im/2 = J* smde id (3.20) 

where d is again the J 1 T = 1 pion-pion scattering phase 
shift [72]. It is convenient to separate this term out explicitly. 
The contribution to the Im J's not included by the two-pion 
intermediate state can be formally added in by writing 

Im J lj2 = JJ >2 sin de id +Im J' li2 . (3.21) 

For example, as can easily be seen by inspection of Fig. III.5(b), 
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the two-pion state does not include the perturbation contri 
bution (order g*, Kg. III.5(c)) to Im J\ it will therefore be 
contained in Im </'. 





FIG. III.5. (a) shows symbolically the amplitudes J and J" 2 ; (b) shows 
the two-pion contribution to Im J lf and Im J 2 , and (c) shows the 
perturbation theory Feynman graph contributing to J^ and J 2 . 

7. Now the development may be completed as follows. Recall 
that the pion form factor satisfied a very similar relationship. 



Im F, = F^ sin 8e id , 



(3.22) 



within the framework of the approximations used in dealing 
with it. If both F^ and J 1? say, satisfy dispersion relations, so 
does [73] /i/j^V 

Now Im J( 2 receives contributions from diagrams similar 
to Fig. III.5(b), but with the two-pion intermediate state 
replaced by any other allowed state. The lightest of these is 
four pions, with a threshold at 16m^, so that Im J' l}2 does not 
vanish above q* = 16m^. In addition, it can be shown from 
the Mandelstam representation that Im J{ 2 is non-zero for 
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). In the region 
Im J' lj2 is determined solely by the perturbation diagram, Fig. 
III.5(c). These results are derived in Appendix B. 

In the region 4m^g 2 ^16m^, then, J* 2 sin de ld is itself real. 
If it is assumed to be real for all g 2 , then 



T 
lm 

Thus one obtains 



^f\ - 
Fj ~ 



The integration here, of course, runs over the region of q r 2 for 
which Im J"{ j2 (#' 2 ) does not vanish. As yet this expression is 
still too e exact' to be useful, since, Im J[ }2 contains everything. 
A not too far fetched approximation now suggests itself. 
Namely one might try to retain in ImJ' l2i (q f2 ) only the 
contribution of order gr 2 , and replace Be F v (q' 2 ) by its lowest 
order value, namely e. The best argument for the validity of 
this approximation is that attempts have been made to 
include more details of Im J^ 2 and Re F v and these refinements 
do not produce significant changes, at least for J 2 . With this 
approximation, then, plus the observation that 



- 



(the subscript BA means Born Approximation), the final 
answer is obtained: 



Inserting this back in the expressions for the imaginary parts 
of the nucleon form factors yields, finally: 

Im JJ >2 ( 2 ) = i \F&*)\* Im J?1 2 (<Z 2 )BA. (3-23) 



The Born Approximation results for Im F\ and Im Fl were 
derived much earlier (see equations (3.16) [74]. 
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It remains, then, to see how well this approximation agrees 
with experiment. The results obtained depend on the P-wave 
TTTT phase shift through F v9 and nothing is yet known experi 
mentally about this phase shift. One might start therefore, by 
asking what kind of a phase shift is necessary to fit the 
experimental results. Prazer and Fulco [71] have used a two 




FIG. III. 6. Pion form factors resulting from various choices of the 

TTTT resonance in the T = U I state. (The graph is taken from 

Frazer and Fulco, [71].) 

parameter resonance formula for the pion-pion scattering, and 
have determined values for the position and width of the 
resonance which give agreement with the experimental 
isotopic vector moment structure as well as with the static 
moment value since they use the unsubtracted form for the 
calculation of Fl(q 2 ) in equation (3.16). Their results are 
shown in Figs. IIL6, 7 and 8. They find that if <5 goes through 
90 at a value q' 2 of the order of 12m^ to 16mf , in a fairly sharp 
resonance, agreement can be achieved. It is significant that 
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to 




FIG. III. 7. A comparison of the proton isotopic vector form factors 
(assuming F^IF^ (0) = F*fF*, (0)) resulting from various choices of 
the TTTT resonance with experiment. The comparison is probably 
meaningful only for small values of g 2 . (Taken from Frazer and 

Fulco, [71].) 



200- 




771-7 



FIG. HI. 8, Two choices for the T = 1 J ITTTT cross-section used in 
fitting the nucleon structure. (Taken from Frazer and Fulco, [71].) 



CALCULATIONS OF NUCLEON FOBM FACTORS 85 

they are unable to obtain agreement without such a resonance; 
its existence, then, may be said to be a prediction of the 
dispersions approach within the approximations which have 
been outlined above. 

It is perhaps even more significant that if a subtracted 
dispersion relation is used for the moment, which should tend 
to improve the validity of the two-pion approximation, a 
resonance is still necessary; in fact, even if two subtractions 
are made a fit to the experiments is impossible without 
appreciable TTTT scattering [75]. 

The basic question is then reduced to whether or not 
experimental or theoretical information about <5 implies the 
presence of this resonance. In retrospect, what has been 
achieved by the analysis of the form factors is a relation 
between two experiments, the nucleon structure and TTTT 
scattering. If eventually TTTT scattering experiments are not 
consistent with the nucleon structure, one will have to decide 
between accepting the assumptions involved in the present 
calculation, and dropping quantum electrodynamics or meson 
theory entirely. 

At present, experimental information on TTTT scattering is 
non-existent though some experiments have been proposed 
[76], Various theoretical calculations have been made or are in 
progress [77]. As an example of the kind of agreement which 
can be obtained through some of these theoretical estimates of 
7777 scattering, together with the approximations outlined 
above, Figs. III. 9 and III. 10 are shown. The theoretical 
curves here are obtained on the basis of equations (3.23) 
combined with <5 obtained from a determinantal calculation, 



assuming a point four-pion interaction of the form 

il 

The agreement seems nice, but is deceptive as the curves have 
been normalized to one at q 2 = 0. The actual zero value may 
be expressed by saying that the theoretically predicted 
isotopic vector anomalous magnetic moment is 0.8 magnetons, 
instead of the experimental 1-85 magnetons. 

Note that the neglect of heavy mass intermediate states, 
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which affect the Im F%s only at large q' 2 , are more likely to 
influence the value of Fl(0) than the shape of Fl(q 2 ). This may 
be seen by returning to the subtracted dispersion relation (3.6), 
namely 



The integrand is now weighted by an extra power of 




-10 - 



FIG. III.9. Pion form factors obtained from an assumed pion point 
interaction A(<p .<p) 2 . (Taken from Baker and Zachariasen, [77].) 

against large g' 2 , but Fl(0) cannot be predicted, and must be 
taken from experiment: 

J?l(0) = 1-8 



An optimist would be satisfied with the prediction of the 
shape, and say that of course inclusion of heavy states might 
allow the evaluation of the moment as well by going back to an 
unsubtracted relation. A pessimist will, on the other hand, 
insist that it is a shortcoming of the present calculation that 
the predicted moment misses by a factor of 2. 

With regard to the charge isovector form factor theoretical 
predictions are again made from equations (3.23). However, 
because of the lack of symmetry between the experimental 
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charge structure of the neutron and proton, it is evident that 
there are sizable scalar contributions to the observed charge 
structures. Therefore experimental information on I{ is 
largely lacking, so comparison with the theoretical predictions 
is difficult. The best information on the neutron charge 
structure, as discussed in Chapter II, is that the charge radius 




-? 2 x10 26 cm- 2 

Fia. III. 10. The predictions of the A(cp .<p) 2 TT-TT coupling for the 
nucleon isotopic vector form factors, superimposed on the experi 
mental results from Fig. II. 5. (Taken from Baker and Zachariasen, 

[77].) 



is ~zero, indicating that for small q 2 the proton charge 
structure coincides with the isovector charge structure. This 
is in agreement with the results of Frazer and Fulco [71]. 

This completes all there is to say at present about the 
isotopic vector form factors. To add to this a discussion of the 
theoretical results which exist on the scalar form factors is 
simple owing to their scarcity. If the unitarity argument is 
accepted, it seems proper to rule out the nucleon-antinucleon 
intermediate state as providing much of the scalar form 
factors. If this philosophy is extended to exclude all heavy 
states, one is presumably left with the three-pion state as the 
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culprit. It is so difficult to even estimate the size or properties 
of what this state gives that the scalar form factor problem 
remains totally un-understood. 

In order that this state explain the observed properties, it 
must provide sizable charge effects but almost no contribution 
to the moments. Some crude estimates of the three-pion state 
contributions, based on perturbation theory [78] and on the 
static meson theory [79] exist, but none have this feature. 
Generally, a scalar static moment emerges which is of the 
order of 50 times too large. 

Moreover, it is disturbing to have to call upon a three-pion 
intermediate state to play a major role in the charge structure 
analysis. One wonders then whether it is a valid approxi 
mation to neglect, say, the four-pion state as has been done on 
the basis of the threshold arguments given in the above 
dispersion approach to the isovector form factors. 

To conclude, then, theoretical understanding of the form 
factors on the basis of meson field theory is fairly meagre; 
totally so for the scalar form factors, less so for the vector ones. 
Predictions, on the basis of certain approximations, can be 
made for the vector form factors; these then agree or disagree 
with experiment depending on whether or not there is an 
appropriate resonance in the TTTT scattering. 



CHAPTER IV 

VALIDITY OF QUANTUM 

ELECTRODYNAMICS AND EFFECT ON 

NUCLEON STRUCTURE INTERPRETATION 

UP to this point the entire discussion has proceeded on the 
tacit assumption that electron-proton scattering, and the 
other experiments mentioned in Chapter II, all are explicable 
in terms of nucleon electromagnetic structure. In Chapter I 
this structure was defined; in Chapter II experimental 
information about it was discussed; in Chapter III theoretical 
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predictions about it were made. As was explained in Chapter 
III, it remains a question whether or not the presently 
fashionable dispersion theoretic calculations actually offer the 
possibility of a quantitative comparison between theory and 
experiment. This question can only be answered by reliable 
calculations of the effects of higher mass intermediate states. 
At any rate the dispersion approach does have the virtue of 
producing a result directly susceptible of experimental test: 
the requirement of a sharp resonance in the J = 1 T = 1 state 
of pion-pion scattering. Evidence of such a resonance is now 
being sought. However even if one is found with the correct 
energy and width to fit F V 2 , the charge radius of the neutron still 
demands a theoretical understanding which at present is 
totally lacking. 

Since there exists then the possibility that pion physics as 
now known [80] may not be able to account for the proton and 
neutron structures inferred from the experimental data, one is 
led to inquire into the possibility of other interpretations of 
these observations. Consider the electron-proton scattering 
experiments; they have yielded the most complete, accurate 
data and their analysis is free of all nuclear physics compli 
cations and pion interaction models. The only alternative to 
the interpretation of electron proton scattering in terms of 
proton structure is in terms of a modification in the present 
theory of quantum electrodynamics of the electron and the 
photon. A finite extension of the electron's charge distribution 
for example would show itself by replacing the proton form 
factors in the Rosenbluth formula (1. 10) by products of an 
electron times a proton form factor [81]. 

There would be two direct and significant consequences of 
accounting for the observed structure wholly or largely by 
electron size: 

(1) In order to achieve agreement with the isovector 
moment distribution it was necessary in (3.19) and (3.23) to 
introduce a fairly sharp TTTT resonance for a total energy of the 
two pions F& 3 -Sm^ in their centre of mass . This was required in 
order to enhance the absorptive amplitude Im jPl(? 2 ) f r small 
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values of q* as seen in Fig. III. 6. This in turn led to a spread 
out distribution, which may be made evident by re-expressing 
(1.11) with (3.5): 



Therefore it follows that, by attributing much if not all of the 
observed structure to the electron, one obviates the require 
ment of the resonant pion-pion interaction. 

(2) The very puzzling difference between the proton and 
neutron charge sizes would disappear. 

What direct evidence is there then for finite electron size or 
for modifications of any sort in the theory of quantum 
electrodynamics? So far, none; indeed it is supported by very 
precise measurements [82], The magnetic moment of the 
electron including radiative corrections is accurate to better 
than one part in 10~ 5 : 

&*pti = 2(1-001162 4xlO- 6 ) 

<7theor = 2(l+ ^ -0-328^) = 2(1-0011596) (4.1) 

and calculations and observations on the Lamb shift are in 
agreement with an accuracy of ^0-2 megacycles [83]. 

These however are low energy experiments and probe 
possible electron size or quantum electrodynamic modifi 
cations only with a very crude yardstick whose unit of measure 
is the electron compton wavelength: H/mc & 3-8 xlO"" 11 cm. 
In order to test the validity of the theory down to distances of 
less than 10~ 13 cm, such as enter into the nucleon structure 
analysis, large momentum transfer experiments are required. 
They contain a small invariant length (corresponding to H 
divided by the invariant 4-momentum transfer) against which 
to measure the small distance behaviour. In fact, at present, 
it is the electron-proton scattering experiments themselves 
which supply the most severe test of quantum electrodynamics 
at small distances since they show that modifications are 
limited to distances of < x 10~" ia cm (corresponding to a root 
mean square radius of 0-8 x 10~~ 13 cm which is conventionally 
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quoted as the proton 'size 3 on the assumption of a point 
electron) [84]. 

How then is it possible to extend the limits of one's know 
ledge of quantum electrodynamics and to decide unambigously 
if the c size' deduced from electron-proton scattering experi 
ments is in fact the structure of the proton or if it is an 
indication of deviations from the theory of photons and 
electrons as now formulated? 

One possibility mentioned earlier is that of using the 
electron proton scattering itself as a consistency check on the 
theory by performing three measurements at different energies 
and angles corresponding to the same # 2 . As remarked by 
Feynman [85] it would no longer necessarily be the case that 
only two functions, F^q*) and F%(q z ) summarize the electron- 
proton scattering data if there are contributions to their 
interaction corresponding to the exchange of other than 
vector photons, or if electrons and protons can virtually 
combine to form a boson. To this may be added the observa 
tion that the same conclusion would also follow if the electron 
were to develop an anomalous moment structure for large g 2 . 

The Bosenbluth formula may be generalized to allow for 
such a possibility in which case (1.10) becomes 



cos 



X 



{[A 



-2/ 1 - - a 2m/ 2 ] [>?- 



sn 



tan* 

where/! and/ 2 are form factors for the electron and m has been 
put equal to zero everywhere except for the one large term 
indicated [86]. Since this formula contains three different 
angular terms, for 2m/ 2 = 0, it allows a more general analysis 
of the observations in terms of three invariant functions of q 2 . 
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It is therefore of great interest to establish whether or not 
the scattering data can be completely summarized in terms of 
only two functions by constructing three (or more) experi 
mental ratios of the scattering at fixed g 2 . If not, one may 
attribute it to an anomalous electron-proton interaction such 
as discussed by Feynman or to an anomalous electron moment 
/ 2 as indicated above. This latter possibility is unlikely to 
show up in the scattering experiments however, because the 
measurements of the electron g value, (4.1), agree to better 
than one part in 10 5 with electrodynamic predictions. There 
fore, unless the anomalous moment / 2 peaks up for large space 
like g 2 , it follows that/ 2 < 10" 5 e/2m and (2m/ 2 ) V/4m 2 <10~ 4 e 2 
for 500 MeV electron energies. 

In order to test for electron form factors or for other possible 
types of breakdown in quantum electrodynamics it is necessary 
to turn to different experiments. It is conceivable for example 
that the photon propagator has to be modified at small 
distances, indicating a departure from Coulomb's law. This 
would mean a change in the Feynman propagator for a photon 
of four momentum q* [87]: 



Another possibility is to alter the electron-positron propagator 
at small distances, or large momenta, in a similar fashion; 
there must in fact be a connection between electron size, or 
form factors, and alterations in the behaviour of the propa 
gator of the electron positron field in order to maintain a 
differential current conservation law [88]. 

Such modifications may be studied directly by electron- 
electron, electron-positron or by photon-electron (compton) 
scattering. The Feynman diagram, Fig. IV. 1, for electron- 
electron scattering shows that it probes electron size and the 
photon propagator as a function of the momentum transfer. 
The diagrams for compton scattering, Fig, IV.2, show that the 
virtual particle in this process is the electron, so that what is 
probed here is the electron propagator at small distances. 
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For these processes to probe the structure of quantum 
electrodynamics at distances <10~ 13 cm it is necessary to 
achieve momentum transfers of ^200 MeV/c. In experimental 
arrangements scattering from target electrons most of the 




FIG. I V.I. Diagram representing the single photon exchange 
contribution to electron electron scattering. 



V 





Fia. IV.2. Feynman graphs for compton scattering: 



energy is in the centre-of-mass motion and the momentum 
transfer is limited by {2mE Q }~ . Therefore a 40 BeV accelerator 
is required in order to realize the necessary momentum 
transfer [89]. As such energies are not now available in the 
laboratory one has two choices either to increase the energy 
available in the centre of mass by replacing the target with a 
colliding beam, or to 'embed' the interesting Feynman graphs 
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of Fig. IV. 1,2 in more complicated processes which avoid the 
centre-of-mass problem by using heavy target particles. 

The first of these alternatives is now being pursued at 
Stanford for electron-electron scattering [90]. Two beams of 
500 MeV electrons are formed in intersecting storage rings and 
their collisions studied by counting pairs of electrons emerging 
back to back in coincidence at various scattering angles. A 
detailed analysis of the radiative corrections to the Moller 
cross-section has been carried out by Tsai [91] for experimental 
parameters of interest here. If there is any disagreement 
between theory and experiment, it will have to be attributed 
to finite electron size or to modifications of the photon 
propagator. Since momentum transfers of up to 700 MeV/c 
will be realized, this experiment, if successful to a 10% 
accuracy, will be an enormously high resolution probe to 
distances of 10~ 14 cm. At present this project is in the 
construction stage and it is hoped that data taking will start 
during 1961. 

The second alternative of achieving large 4-momentum 
transfers at present laboratory energies by using a heavy 
target particle is also being pursued [82], [92]. In particular, 
RIchter has measured large angle production of electron- 
positron pairs by photons incident on hydrogen; further 
experiments aiming for greater precision and higher energies 
are now in progress. These serve as probes of quantum 
electrodynamics in the following way. If quantum electro 
dynamics is valid it is possible to summarize the proton's 
structure in terms of measured F^q*) and F%(q 2 ). This is seen 
in the Feymnan diagrams for pair production, Fig. IV. 3. 
These graphs show that only the electromagnetic vertex of a 
physical proton appears. Therefore any disagreement between 
theory and experiment can be interpreted only in terms of 
modifications in quantum electrodynamics, and a comparison 
between theory and experiment yields new information on the 
small distance behaviour of the virtual intermediate electron- 
positronline. The present agreement between theory and experi 
ment indicates no modification in the electron propagator at 
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distances ^10~ 13 cm. Further experimental refinements are 
anticipated to lower this limit by a factor ^3. 

These large angle pair studies provide direct information 
only on the electron propagator. The photon propagator and 





FIG. IV. 3. Feynman graphs for pair production: 



electron structure appear together with the proton structure 
as in the elastic scattering [93]. 

In order to investigate electron size and the photon 




. IV.4. A Feynman graph, contributing to trident production: 



propagator further consider another experimental arrangement 
in which the pair is produced by an incident electron, i.e. 
trident production. Fig. IV,4 shows a typical Feynman graph. 
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By tying the photon of Fig. IV. 3 on to an incident electron 
the photon propagator and electron vertex are introduced. It 
has been shown [94] that the ratio of this to pair production 
provides a sensitive test of these structures and hence of the 
interpretation of the electron-proton scattering experiments 
in terms of proton structure. 

In conclusion then three points may be repeated: 

(1) It remains for experiment to verify unambiguously the 
analyses of the electron-proton scattering and of the other 
experiments discussed in terms of nucleon structure in 
Chapter II. 

(2) It remains for theory to verify that the theoretical 
approximation scheme is appropriate for a quantitative 
comparison with experiment. 

(3) If (1) and (2) are fulfilled, it remains finally for 
experiment to verify the TTTT resonance in the J 1, T = 1 
state predicted by the dispersion analysis of the isovector form 
factors. 

APPENDIX A 

It is desired to verify the relation used in equation (3.7) or 
more generally, to prove the identity 




-i ^x(^\(^(x}) + \^)(-i7J x -m) -L-e-**-*. (A.I) 



The notation here needs some explanation: |ot (+) ) and (^ ( ~ } \ 
are arbitrary Heisenberg states with the indicated boundary 
conditions; |aj9 (+) > is a state containing, in addition to 
whatever is contained in |a (+) ), a nucleon of momentum p. 
is an arbitrary Heisenberg operator, and ip(x) is the renor- 
malized Heisenberg nucleon field operator. In the first term, 
the symbol d^ np is meant to be zero unless the state [/?} 
contains the nucleon of momentum p, and if |/?) does contain 
the nucleon p, the state \n) is to represent everything in |/S) 
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other than this nucleon. The first term in equation (A.I), 
then, represents the possibility of the nucleon p not interacting 
with anything, but just passing from the initial to the final 
state as a free particle. 

The proof of equation (A.I), which is generally referred to 
either as the 'contraction rule/ or the 'reduction formula/ is 
not difficult, and is sketched in what follows [95]. 

Consider the matrix element lira <j8 ( " } |(Oy(a;)) + |a (+) >. Then 

lim <<~ ) 



= 2X/3 ( - } |0|w<+>> lim <<+>|y(3!)|a<+>). (A.2) 

n x<r^ 

Now, since ip(x) -+ ^ in (#) as X Q -+ oo, 



lim 

Q-+ 00 



It should be noted that no factors of Z 2 appear here since ip is 
the renormalized field. If (A.3) is substituted into (A.2), and 
the Fourier transform is inverted, one obtains 



= lim 

X -+- 

Now observe the identity 

f(x) = -i^yf(y}(iy y ~-m)K+(y, x), 

which holds since the nucleon propagator K + satisfies 
(iK/ y m)K + (y, x) = id*(yx). 



A partial integration may be performed in this identity, 
yielding 

f(x) = -id*yf(y)(-W v ->)Z+(y 9 x) 

^^ (A. 5) 
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Equation (A,5) may now be applied to 
equation (A. 4). Using the fact that 



n 



u m 



-ip.y 



one finds 




. (A.6) 



The first term in (A.6) may be rewritten as sketched below. 
Write 
lim <^ < - ) |(Oyi(2/)) + |a'+)> 

= lim J dS'-'lvK^K-' ><> ( -> 1 0|a<+)}; 

2/ ~>-+co w 

now lira ip(y) = y5 out (^), so that the same algebraic manipu- 

y -H. + CO 

lations as were used to arrive at (A. 4) from (A. 2) may be 
used to replace the first term in (A.6) by 

2<5 KIJ /n<->|0|a<+>>. 

n 

This completes the verification of equation (A.I). 
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In this appendix, a more complete discussion of the regions 
in which Im J' is non-zero will be presented. This is equivalent 
to discussing the singularities of the annihilation amplitude for 
the reaction TT +TT ~> N +N. It will facilitate the discussion to 
assume, as has been done before, that nucleons have spin zero. 
The amplitude for this reaction is then a single scalar function 
of two variables, which may conveniently be chosen to be 
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Here the two meson 4-momenta are denoted q l and g 2 ; those 
of the nucleons p and p z . Thus the amplitude for the anni 
hilation may be written 



In terms of s and t, the angle of scattering in the centre-of-mass 
system is 



= 
C S 



A particular partial wave amplitude for the process is then 
obtained by l 

A^s) = \ A(s, *)P z (cos 6) c(cos 0); (B.3) 



the integration over d(cos 0) may be replaced by one over dt by 
virtue of equation (B.2). It is one of these partial wave 
amplitudes, specifically the 1 = 1 amplitude, which is the 
analogue of the functions J and / 2 discussed in Chapter III. 
The fact that the nucleon spin has been ignored is responsible 
for the appearance of only one function here instead of two; 
however the kinematics have not been affected, so that the 
singularities of A 1 which will be obtained here are identical 
with those of the actual functions J t and 7 2 . 

If the contraction rule is applied to the annihilation 
amplitude, then as in the derivation of equation (3.8), one 
obtains 



Equation (B.4) may be thought of as the decomposition of the 
TT+TT -+N+N amplitude into the processes TT+TT -> n -+ 
N+ff (the first term) and 77 -> n+N followed by ir+n -+ N 
(the second term). Now, the first term will be singular if 

PU 

that is, if 



100 



APPENDIX B 



Since the lightest state n here consistent with the various 
conservation laws is two pions, the first term is singular when 



s = 



The next heaviest state is four pions, which produces a 
singularity for 



thus for 4m^s<16ra^ all of the singularities in the first term 
of (B.4) are due to the two-pion intermediate state, pictured in 
Kg. III.5(b). 

The second term in (B.4), on the other hand, is singular if 



in other words, if 



PI = fe-^i) 2 = t. 



Here the lightest allowed state is one nucleon, for which 
P| M 2 ? an d which therefore produces a pole at t = M*. In 
diagrams, this is represented by Kg. III.5(c). The next state 
is one nucleon and one pion, for which P^(M +mf) 2 , so there 
are singularities from the second term of (B.4) for 



This term is illustrated in Fig. (B.I). 




FIG. B.I. Graph illustrating the one nucleon one meson intermediate 
state contributing to the second term of equation (B.4). 
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Now singularities will occur in A^s), according to equation 
(B.3), either when A(s, t) is singular as a function of s, or when 
the integration over cos 6, or equivalently t, includes a region 
of t for which A(s, t) is singular in t. It is the assumption of the 
Mandelstam representation that no further singularities exist 
in A AS). 

A AS) will then have singularities for s>4m^, of which those 
between 4m* and 16m^ are due only to the two-pion state 
contribution to the first term of (B.4). It will also have 
singularities whenever the range 

_i 

is consistent with t = M 2 or t^(M +mj 2 . It is possible to 
have t M 2 in this range when 



it is possible to have t^(M 4-mJ 2 when 



The singularities of A AS) between and 4m^(l ra^/4.34" 2 ) are 
therefore all due to the single nucleon pole in the second term 
of (B.4). 

Finally, as has been indicated in Chapter III, Im^l z ($) is 
proportional to the discontinuity of A AS); thus ImA^s), or 
equivalently Im J lj2 of the actual problem will be non-zero if 
^4ra^ 3 or 5^4m^(l -m^/4Jf 2 ). Since Im J(^ is the same as 
Im J 12 with the two-pion state omitted, Im J{^ will be non 
zero for 5>16m or 



NOTES 

[1] The following notation will be consistently used: 4- vectors are 
denoted a, or a^ or (a , a). The scalar product is defined as 

a . b = &O&Q ** ^ == a /$ti> 

raising and lowering of indices will be suppressed. The Dirac 
matrices are defined by 

y^ = ()8, j 
and 
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The symbol @ will be used for a . y = a^y^. The operator V ; , is 

/a \ 

I ^7 , V I ; thus the free Dirac equation is written (H = c = 1) 
\oi } 

(iW-M)ip 0. 

The spinors which are used are normalized to u^u v = 2M for 
positive energy solutions where u == u*y Q , and VJQ^ = 2M 
for negative energy solutions. All field operators which appear 
are renormalized Heisenberg fields; Heisenberg states are 
denoted \p), for example, for a single particle of momentum p. 
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[3] S. Schweber, H. Bethe, and F. de Hoffmann, Mesons and Fields 

Vol. 1 (Bow Peterson and Co., Evanston, 111. 1955), p. 119. 
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[6] The experimental basis for these definitions is this. The zero 
energy limit of the amplitude for Compton scattering, i.e., the 



Thomson limit, is proportional to V_ } e z /M. This comes 
from the two Feynman diagrams 






which for zero photon energy measure just the product of two 
vertices connecting physical (because lim (p&) 2 = M 2 ) 



proton states and also measure ^(0) since q 2 = for real 
photons. jP 2 (0) is determined by measurement of the proton 
magnetic moment in an external magnetic field. 

[7] J. Hamilton, The Theory of Elementary Particles' (Oxford 
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(r\ iN = 8J7f (g-'M^o/JFf (0) 
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e4<f> V /{pjylJjylnXttlJglflp) 

vW " n \$n-E 9f -E n -K 



_ 

- 

where j N (iyjMty. The pole for the one neutron inter 
mediate state in Fig. 11.10 arises from the second term for \n)a 
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NOTES ADDED IN PROOF 

1. Page 66. R. Cutkosky (Tenth Annual International 
Conference on High Energy Physics at Rochester) has 
shown that such a replacement generally produces the 
absorptive part of any Feynman graph. 
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2. Page 85. Frautschi and Walecka [Phys. Rev. (1960) in 
press] and Bowcock, Cottingharn and Lurie [Nuovo 
Cimento, 16, 918 (I960)] have included the effect of TTTT 
scattering on irN scattering, and find no significant effect 




on the large 33 phase shift. However Frautschi [Phys. 
Rev. Letters, 5, 159 (I960)] has shown that with the 
Frazer-Fulco [71] parameters for the TTTT interaction, un 
acceptable small P-wave phase shifts are found in the irN 
problem. Bowcock et al., on the other hand, determine 
the TTTT interaction by fitting the energy dependence 
of the S-wave nN phase shifts and the radius, <r 2 ) 2 , of the 
magnetic form factor F\, and find the TTTT resonance at a 
higher energy (around 5m J and narrower in width than 
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that of Frazer and Fulco, leading to too large a static 
magnetic moment by a factor of the order of two. These 
results indicate that the two-pion approximation may 
well be inadequate. 

As yet there is no clear experimental evidence for or 
against the existence of the -mr resonance. 
3. Page 19 and page 87. The latest experimental informa 
tion on the proton form factors, as obtained from electron- 
proton scattering, is shown in the graph taken from 
Hofstadter, Bumiller and Croissiaux. These curves are 
in agreement with data by Beckelman, Cassels and 
Wilson, and by L. Hand. (Report to Tenth Annual 
International Conference on High Energy Physics at 
Rochester.) 

It is interesting that FI appears to be levelling off, 
indicating that approximately 40 per cent of the time a 
proton is a "proton." It is also worth noting that the 
shape of -F 2J together with the results on page 64, shows 
that Im jP 2 cannot be positive definite. 

L. Hand, Phys. Rev. Letts., 5, 168 (1960). 

Bumiller, Croissiaux, and Hofstadter, Phys. Rev. Letts., 5, 261, 

263 (1960). 
Wilson, Berkelman, Cassels, and Olson, Nature, 188, 94 (1960). 




112515 



